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ABSTRACT 

NGC 891 is a nearby edge-on galaxy that is similar to the Milky Way and has a hot X-ray emitting 
halo that could arise from accretion, a galactic fountain, or a combination of the two. The metallicity 
of the gas can help distinguish between these models, and here we report on results that use 138 ks 
of archival Chandra data and 92 ks of new XMM-Newton data to measure the temperature and 
metallicity of the hot halo of the galaxy. We find good fits for a thermal model with kT ~ 0.2 keV 
and Z ~ O.IZq, and rule out solar metallicity to more than 99% confidence. This result suggests 
accretion from the intergalactic medium as the origin for the hot halo. However, it is also possible to 
fit a two-temperature thermal model with solar metallicity where kT\ ~ 0.1 keV and kTi ~ 0.25 keV. 
A consideration of the cooling rate and scale height prefers the single-temperature model. We also 
find that the cooling rate in the hot gas cannot explain the massive H I halo in the steady state. In 
addition, a galactic fountain model cannot eject enough mass to account for the H I halo, and we 
speculate that the neutral halo may be gas from a prior outflow that has since cooled. 

Subject headings: galaxies: individual (NGC 891) - galaxies: abundances - galaxies: halos - X-rays: 
galaxies 



1. INTRODUCTION 

The hot halos of spiral galaxies like our own are reser- 
voirs of 1 — 3 x 10 6 K gas whose origin is unclear. These 
halos are potentially important to galaxy formation the- 
ories because some closed-box models of Galactic chem- 
ical evolution predict a larger number of metal-poor 
stars in the solar neighborhood t han we observe (the 
G-dwarf or M-dwar f problem, see Ivan den Berghlll96l 
iWoolf fe~W cst 2012 ), and this issue i s not isolated to the 
Milky Way (e.g. iWorthev et aLl ll996). This problem can 
be solved by the accretion of low -metallicity gas onto the 
disk from the halo. In addition, Sancisi ct al. (20QH) ar- 
gue that about lM yr _1 of accretion onto the disk is 
required to sustain the star form ation rate (SFR) of typ- 
ical spirals; iBinnev et al.l ([2000) find that the SFR in the 
solar neighborhood has been nearly constant for the past 
HGyr for a Salpeter initial mass function, so star for- 
mation may indeed be sustained for long times. It is not 
yet clear whether the required accretion occurs primarily 
via condensation of a radiativcly cooling hot halo or cold 
flows and minor merge rs (for a review of cold accretion, 
see lSancisi et alJl2008f ). 

It is also unclear whether the hot halo is gas accreted 
fro m the local intergala ctic medium (IGM) as suggested 
bv ICox fc Smith! (|1976f ). Instead, they may be gas ex- 
pelled fro m the disk by sup ernov ae (SNe) in a "g alactic 
fountain" ([Shapiro fc Field! [l97l iBregmad 11980 ). Both 
can explain the observed hot halos: in the case of accre- 
tion, infalling gas from the IGM shocks and is heated to 
the dynamical temperature, whereas a galactic fountain 
forms the hot halo as superbubbles powered by hundreds 
of SNe break out of the galactic disk and vent their con- 
tents to the halo. 

Since both the galactic fountain and accretion from 
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the IGM are expected to occur, it is not clear which con- 
tributes most to the hot halo. For example, if the super- 
nova rate is too low, the bubbles may never break out 
or only contribute to a halo with a small scale height. 
As both mechanisms can produce a halo with a simi- 
lar temperature and disk condensation rate, the simplest 
way to distinguish them is the metallicity of the hot gas, 
which should be close to the IGM value in the accretion 
scenario (~ 0.05 — 0.3 . 2^ in a spiral-rich cluster such as 
HCG 16: iBeTsole et all 120031 and solar or super-solar in 
the fountain scenario. 

The hot halos are visible in the soft X-rays where the 
metallicity can be measured directly. Since edge-on spi- 
rals allow for clean separation of the disk and halo emis- 
sion while amplifying the surface brightness of the halo, 
they are the ideal candidates for this measurement. 

In this paper, we present an X-ray study of the hot 
halo of NGC 891 using Chandra and XMM-Newton data 
to measure t he halo metallicity . NGC 891 is a nearby 
(d ~ 10 Mpc; I Temple et al.l l2005D edge-on spiral thought 
to be a close M ilky Way analog in luminosity and spec- 
tral type (e.g. Ivan der Kruit fe Searld E981), and, as 
in our Galaxy, there is dynam ical evidence for a bar 
(IGarcia-Burillo fc Guelinl fl995h . NGC 891 is not inter- 
acting with other galaxies, is almost p erfectly edge- 
on, and has no o ptical or H I warp (|Rupenl 119911 : 
ISwaters et al.l|j"997l ). allowing an unambiguous distinc- 
tion between halo and disk em ission. However, it i s 
known to be slightly lopsided (jSancisi fc Allen] fl979h . 
The primary aim of our study is to measure the metal- 
licity of the hot halo directly from the X-ray data. A 
summary of basic parameters for NGC 891 is given in 
Table [TJ 

NGC 891 is of particular interest be cause it is the X-ray 
brigh test of normal edge-on galaxies (Bregma n fc Houckl 
1997) and has been held up as a prime example of a 
galactic fountain in action. The evidence in favor of 
such a fountain is quite strong, enabling us to deter- 
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Fig. 1. — Left: Composite 2MASS H-band image of NGC 891 (green) with bright emission from the adaptively smoothed 0.3 — 2.0kcV 
Chandra image (red). Right: Same as left, but with 0.3 — 1.2 keV combined adaptively smoothed XMM-Newton image in red instead. Te 
dark holes are point source masks shown in Figure [2] (the large one is the ULX) . 




Fig. 2. — At left we show the combined EPIC 0.3 — 1.2 keV image with point sources removed before exposure correction. The remaining 
three panels show the Swiss-cheese point-source and chip-gap/bad pixel masks created with XMM-ESAS, centered on the galaxy. The 
XMM-Newton image in Figure [l] is a smoothed, exposure corrected version of the counts image shown here. 



mine whether the X-ray emitting gas is an important 
participant in the fountain or largely decoupled accreted 
material that slowly settles onto the disk. We review the 
evidence presently. 

First, the disk appears to be dynamically connected 
to the halo. The cold H I halo rotates differentially, 
with rotation velocity decr easing with distance from the 
disk (jOosterloo et al.ll2007l) . and the spatial structure of 
the radio emis sion correlates w ith that of the Ha halo 
(lDettmarlll990D and th e X-rays (jBregman fc PildisHl99l 
Strickland et al.l 1 20041 ). being concentrated in the north 
and near the disk. Th e SFR is also highest here, which 
iTiillmann et~all (f2006l ) argue leads naturally to the ob- 
served distribution since the amount of hot gas is propor- 
tional to the energy injecte d by supernovae. This incon- 
sistent with the dust halo ()Howk fc Savagdll997l 120001) . 
which shows more pronounced filaments on the north- 
ern side (although perhaps just as many on the southern 



side) , but we note that iGarcia-Burillo et al.1 (Il992t ) find 
that the extended molecular halo is consistent with being 
symmetric about the rotation axis. 

Second, there is strong evidence for magnetic fields lac- 
ing the halo that are naturally explained by a galactic 
fountain. There is radio synchrotron emission cospa- 
tial with the H I halo out to lOkpc away from the disk 
(jOosterloo et al.ll2007| ) , and there are collim ated Ha fila- 
ment s that extend a few kpc above the disk (jRossa et al.l 
l2004f h Thi s suggests that the magnetic field originates 
in the disk (|Dahlem et al.lll994l ). 

Finally, cle ar evidence o f supe rnova activity is visible 
in the disk. IRossa et al.l (|2004f ) find superbubbles and 
shells near the midplane, as well as filamentary structures 
that appear to be connected to the disk superbubbles. 
These may signify breakouts from the disk where the 
hot medium can rise into the halo. 

The rest of this paper is organized as follows. In Sec- 
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TABLE 1 

Basic Parameters for NGC 891 



Parameter 


Units 


Value 


a (J2000) a 
A ( T9nnrA a 

O ^JZUUUJ 




02h22m33.4s 




-L-AOA Cirri 


d h 


(Mpc) 


9.8 






0.00176 


i d 


(deg) 


89 


Scale 3 


(kpc arcmin - x ) 


2.2 


Morphology 


SBb 


M m d 


(M ) 


1.2 x 10 9 


M tot d 


(M Q ) 


1.4 x 10 11 


Lx f 


(erg s x ) 


4.4 x 10 39 


mv a 


(mag) 


10.81 


SFRS 


(Moyr- 1 ) 


3.8 



References. — (a) NED (b) 

ITikhonov fc Galazutdinoval J2005I) (c) 

Id^^ aucoule urs et aL| l|19 9l|) (d) | Qpst erloo et all 
120071) (e) IGarcia-Burillo fc Guelinl 119951) (f ) 
IBregman fc Pildij 111994) (g) IPopescu et al.l 120041) 

tion [21 we describe the X-ray observations of NGC 891 
and how we processed the resultant data. Section [3] con- 
tains the temperatures and metallicities we derive from 
fits to the spectra, as well as estimates of the uncertain- 
ties in these parameters from various systematic effects. 
We present two competing models that cannot be spec- 
troscopically distinguished with the present data, so in 
the first part of Section|4]we appeal to the expected cool- 
ing rate and morphology of the hot halo in each model. 
Afterward, we consider more generally whether the hot 
gas is coupled to the neutral halo and whether the sys- 
tem is in a steady state. We conclude by summarizing 
our results and arguments. 

2. X-RAY OBSERVATIONS k PROCESSING 

We processed new and archival X-ray data from the 
Chandra X-ray Observatory and XMM-Newton (Ta- 
ble [T]). Our primary goal was to measure the metal- 
licity of the hot halo, which is visible in the soft X-rays. 
To do so, we extract spectra from regions on the detec- 
tors corresponding to the halo and fit them with a ther- 
mal plasma model, assuming collisional ionization equi- 
librium. In the remainder of this section, we describe 
the observations and processing, surface brightness maps, 
and our spectral fitting method. 

2.1. Chandra Reduction 

We used two archival Chandra observations from 2000- 
11-01 (obsID 794) and 2003-12-10 (obsID 4613) with a 
combined exposure time of 170 ks. These da t a have pre- 
viously been publish ed by iStrickland et al.1 ()2004h and 
iTemple et al.1 (|2005[ ). The data were reprocessed us- 
ing standard Chandra Interactive Analysis of Observa- 
tions (CIAO v4.4) software recipes to create level=2 files. 
We extracted light curves from background regions in 
the 0.3 — 10 keV range to look for flaring exceeding 3cr 
from the mean count rate. A large flare affected the 
last 23 ks of obsID 794, and intermittent short flares 
affect about 5ks of obsID 4613, leaving 138 ks of good 
time. Our analysis is confined to the exposure-corrected 
ACIS-S3 chip except where otherwise mentioned. We 
then used the wavdetect algorithm implemented in 
CIAO to create point source masks in the 0.3 — 8keV 



Fig. 3. — Red: Inner halo from 0.3 — 2.0 keV adaptively smoothed 
Chandra image. Green: 2MASS -ff-band image. Blue: 0.3— 1.2kcV 
adaptively smoothed XMM-Newton image showing the outer halo. 
Contours: Total H I column density, where contours are 10 19 10 20 , 
10 21 , and 10 22 cm - 2 . The Galactic column towards the galaxy is 
approximately 6.5 X 10 20 cm -2 . 



bandpass, with manual filtering for anomalous regions 
near the chip edges. In practice, this eliminates point 
sources down to a flux of F x ~5x 10 -16 ergs -1 cm -2 
in the 0.3 — 8kcV band near the center of the merged 
exposure (where the PSF is small) based on the faintest 
sources detected there. The 5er point source sensitivity 
in the 0.3 — 8keV band for the merged exposures is about 
F x « 3.8 x 10 -16 ergs cm -2 . 

2.2. XMM-Newton Reduction 

We obtained a 133 ks XMM-Newton observation of 
NGC 891 on 2011-08-25 (obsID 0670950101 ). An archival 
shorte r exposure has been published by ITemple et al.l 
(2005), but it is severely affected by flaring and shows 
diffuse emission only near the bulge, so we omit it from 
this study. The data were obtained using the "medium" 
filter, meaning at 1 keV the total EPIC effective area is 
about 1200cm -2 (cf. 555cm -2 for Chandra). The data 
were processed using the standard reduction tasks for the 
EPIC cameras in the Scientific Analysis System software 
(SAS vll.0) and screened for proton flaring, after which 
92 ks of good time remains. Unlike the Chandra images, 
this exposure is severely contaminated by the presence of 
a bright new source with an absorbed 0.3 — lOkeV flux of 
Fx w 1.6 X 10 -12 ergs -1 near the disk of the galaxy. This 
source, which we ar gue is an ultraluminous X-r ay source 
(ULX) in NGC 891 (|Hodges-Kluck et al.l [2012] ) . was not 
present in the Chandra observations nor any archival X- 
ray observation publicly available. 

After standard processing, including exposure correc- 
tion, we created Swiss-cheese point-source masks us- 
ing the cheese and cheese-bands scripts in the XMM- 
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TABLE 2 

X-ray Observations of NGC 891 in this paper 



Date 


Observatory 


Detectors 


Exp. Time 


Sum GTIs 


Obs. ID 










(ks) 


(ks) 


1993-08-27 


ROSAT 


PSPC 


31 


31 


RP500266N00 


2000-11-01 


Chandra 


ACIS-S 


51 


30 


794 


2003-12-10 


Chandra 


ACIS-S 


121 


108 


4613 


2011-08-25 


XMM-Newton 


EPIC MOS/pn 


133 


92 


0670950101 




0.010 



0.002 



0.000 








; East 




0.008 


- □ Chandra 




- A XMM Long 


c 




] O XMM North 


L 




- * XMM South 


o 


0.006 


- CXO South 


o 






\ 

CO 






in 


0.004 




u 






CO 







z (arcmin) 

Fig. 4. — The bottom plot shows surface brightness as a function 
of height z above the galactic midplane in the Chandra 0.3— 2.0keV 
data (black squares) and XMM-Newton 0.3 - 1.2 keV data. The 
adaptively smoothed Chandra image is shown with the extraction 
boxes overlaid. Point sources have been subtracted, but all diffuse 
emission is included, and the values are blindly corrected for the 
point source masks by scaling the true enclosed area to the box 
area. Exponential fits are overlaid. Sec text for discussion. 



Newton Extended Source Analy sis Softwarj^ (X MM- 
ESAS), which are included in SAS (jSnowden et al.ll2004j) . 
We used a cutoff flux of F x — 5 x 10 -16 ergs -1 cm -2 
in the 0.4 - 7.2 keV bandpass. The XMM-Newton 
point source sensitivity for this exposure is about 4 x 
10 -16 ergs -1 cm -2 in the pn and slightly worse in the 



MOS detectorsE 

When compared to the masks produced for the Chan- 
dra data, we find good agreement near the central region 
except for the ULX, which we excise to the 90% encircled 
energy radius (1 arcmin). We also note that chip 5 on the 
MOS 1 detector has anomalously high background and 
is removed from analysis; however, no significant source 
emission falls on this chip. 

2.3. X-ray Images & Scale Height 

The diffuse emission associated with the hot halo can 
be seen in adaptively smoothed maps of the soft emission 
with point sources removed (Figure [TJ the large cavity on 
the right is the ULX mask, and the point source masks for 
each EPIC camera are shown in Figure [5]). The left panel 
shows the bright emission close to the disk seen by Chan- 
dra. Little emission comes from the disk itself (where the 
soft X-rays from the ISM are absorbed by the high Nh 
column), and the emission immediately outside is dom- 
inated by the bulge near the galactic center. The non- 
bulge emission is primarily seen to the north of the galac- 
tic center and is concentrated within about 1.5 kpc of the 
disk. Thi s is the halo emission mo st obvious at lower sen- 
sitivities (Brcgman & Pildis 199 4]), and it cor relates with 
the diffuse ionized (Ha) halo (|Dettmar|[l990l ) and higher 
SFR on the northern side. In the deeper XMM-Newton 
image to the right, it is evident that the diffuse X-ray 
emission extends much farther from the bulge (also de- 
scrib ed in the ROSAT observation in lBregman &; Pildisl 
H99l and is more symmetric on larger scales. The dif- 
fuse emission is much more evident in the XMM-Newton 
image in Figure [T] due to the higher effective area of the 
telescope and the moderate angular resolution (6 arcsec) 
compared to Chandra (0.5 arcsec). 

The difference between the emission on large and small 
scales means that it is possible that the physical prop- 
erties of the emitting gas are also different. For fit- 
ting purposes it is also convenient to extract spectra 
from regions near the disk (where the H I absorption 
is strong and variable, see Figure __]) and regions where 
the dominant Nh column is the Galactic one. Since Nh 
is somewhat degenerate with kT and Z at CCD spec- 
tral resolution, the cleanest measurement of Z can be 
made in the latter case. We therefore divide the X- 
ray emission into two regions: the "inner halo" (where 
absorption must be modeled explicitly) and the "outer 
halo" (where we use the Galactic column only). The in- 
ner halo is defined as the region more than two optical 
sc ale heights above the disk midplane (H ~ 200 — 400 pc 
in lXilouris et al.lll998l and we adopt 300 pc) but below 
where the H I absorbing column falls to half the Galactic 



'• see http: / /heasarc. gsfc.nasa.gov/docs/xmm/xmmhp_xmmesas. html 3 for nominal values, see http: / /xmm.csa.int/external/xmm_user_support/doci 
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value of N H w 6.5 x 10 20 cur 2 (|Kalberla et al.l [2005h . 
which is between 1.5-3 kpc depending on position along 
the disk (Figure [3]). The inner boundary is chosen to 
exclude X-ray binary emission, as ISM disk emission is 
negligible. The outer halo is simply the region above this 
(where the Galactic column dominates absorption) out to 
where the source emission falls below the background. 

In practice, defining the rotation axis of the galaxy 
as z and the projected midplane axis as f, the in- 
ner halo has boundaries 0.3 arcmin < z < 1.3 arcmin 
(0.6kpc< z <2.8kpc) and the outer halo z > 1.3 arcmin 
[z > 2.8 kpc). For extraction purposes, we make boxes 
whose length in f is the projected disk length (7 arcmin 
or ~15 kpc). 

Filamentary structure and signs of outflows are com- 
mon in X-ray images of edge-on (mostly starburst) galax- 
ies (e.g. iStrickland et al J 120041 ). We searched for struc- 
ture by making unsharp masks of the Chandra and 
XMM-Newton data and looking for overlap. Although 
apparent filaments are visible in each dataset, there is 
little agreement between images. On closer inspection, 
the putative filaments often include emission near point 
source masks. Because the PSF wings produce enhanced 
surface brightness over the background while the mask 
itself is a decrement, filamentary structure can be mim- 
icked when a subtracted source is near two (intrinsically 
uncorrelated) regions with higher surface brightness. We 
also find that several of the structures in the XMM- 
Newton map include point sources that were too faint 
or distorted to be identified by the cheese algorithm but 
are obvious in the Chandra data. We find no evidence 
for outflow-like or filamentary structures in the data, but 
as the emission is integrated through the halo we cannot 
rule it out. 

The outer halo is approximately symmetric about z 
and f and is best seen in the XMM-Newton image (Fig- 
ure [1]), as the outer halo is about as wide as the Chandra 
ACIS-S3 chip. To estimate the scale height, we use the 
exposure-corrected, point-source subtracted data in the 
0.3 — 2.0 keV bandpass for Chandra and 0.3 — 1.2 keV for 
XMM-Newton. 

If circular annuli are used, a blind radial profile in- 
cluding the bulge emission in the Chandra image yields 
a surface brightness scale height of ~2 arcmin (~4 kpc). 
The bulge and halo are separate components, with the 
bulge having a smaller scale height and larger surface 
brightness. This suggests a larger scale height for the 
halo. 

As the halo may not be spherically symmetric, we can 
measure the profile in z using regions parallel to the mid- 
plane of the galaxy (simila r to the approach used for the 
NGC891 Chandra data in IStrickland et al.ll200l . This 
profile is shown in Figure 01 The image on top is the 
adaptively smoothed Chandra image, with the 6 arcmin 
by 0.25 arcmin boxes we used overlaid on the image. 
The mean surface brightness in each region is shown in 
the bottom panel, and is blindly corrected for the point 
source masks and other unexposed areas by scaling to 
the nominal area in each box. 

The Chandra data are shown as the black boxes in the 
bottom panel of Figure 01 and includes the bulge emis- 
sion. To derive the scale height, we fit an exponential 
profile (I(z) = Joe + C), finding zo ~ 0.6 arcmin 



(1 .3 kpc) for the Chandra data. This is similar to that 
in IStrickland et al.l (|2004D based on their encircled light 
values, an d somewhat smaller than the ROSAT value re- 
ported bv lBregman fc Pildisl (|1994[ ). The XMM-Newton 
values are shown as red triangles in Figure [4] and result 
in a slightly larger zq ~ 0.9 arcmin (2 kpc). As these val- 
ues are extracted in wider boxes (6 arcmin by 0.5 arcmin) 
because of the large point source masks, we expect a bias 
towards higher Zq. On the other hand, the XMM-Newton 
images mask most of the bulge emission (under the ULX 
point source mask) that is included in the Chandra fit. 
Since the spectrum of the bulge emission shows that it is 
clearly a different physical component (which we expect 
to have a smaller scale height), the Chandra value quoted 
above is probably too small. Since the bulge emission is 
not entirely masked in the XMM-Newton image, we can 
improve this value by ignoring a region within about an 
arcminute of the galactic center. We do this by extract- 
ing profiles in 3 arcmin by 0.5 arcmin north and south 
of the bulge in the XMM-Newton data, seen as green 
diamonds and blue stars in Figure [4] respectively. 

The similarity between the values for the long boxes 
and the northern shorter boxes in Figure 2] shows that 
the correction is small, with the dominant component 
being the bright emission near the disk on the north 
side. The halo emission itself is asymmetric both north- 
to-south and east-to-west. The north-south asymmetry 
is well known and can be easily seen in Figure HI where 
the boxes are centered on the galactic center. The east- 
west asymmetry can be seen in the horizontal dashed 
and dotted lines in Figure 01 which represent a constant 
"background" in the fit to the Chandra data. Despite 
nearly identical scale heights on either side of the galaxy, 
the west side is decidedly brighter out to the extent of 
the S3 chip. Emission in the wings of the PSF of bright 
subtracted point sources can account for less than 20% of 
the discrepancy. A cluster of faint point sources cannot 
be ruled out. 

On the north side of the bulge, we find z$ ~ 1 arcmin 
(2 kpc) east of the galaxy and zo ~ 0.6 arcmin (1.3 kpc) 
on the west side. In contrast, Z the surface brightnesses 
obtained on the south side of the bulge indicate scale 
heights of zo ~ 1.5 arcmin (3.2 kpc) east of the disk and 
zo ~ 1.8 arcmin (3.9 kpc) to the west. The difference 
in these values suggests that the inner halo may be a 
separate physical component. 

However, the scale height is very sensitive to the in- 
nermost surface brightness values. These values are the 
most affected by the point source masks (which are more 
concentrated near the disk), the bulge, and the residual 
10% of the ULX emission. The change in surface bright- 
ness across each extraction zone is also highest near the 
disk, so our blind correction to the nominal area increases 
the uncertainty in these bins far beyond the y/~N errors 
plotted in Figure 2) 

Still, it is likely that the true scale height, at least in 
the outer halo, is even larger than the Zq ~ 1.5 arcmin 
found in the southern boxes in the XMM-Newton data. 
When we correct for the bulge and ULX by subtracting 
model profiles (assuming a circular PSF for the ULX and 
a circular model fit to the bulge in the Chandra data) , zq 
increases to near 1.0 arcmin in the northern boxes (where 
the inner halo is bright) and 2.5 arcmin in the southern 
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0.3 1.0 3.0 

Energy (keV) 

Fig. 5. — Chandra inner (red) and outer (black) halo spectra with 
model components for the inner halo spectrum overplotted. The 
data have been rebinned to 5<r S/N per bin for plotting purposes. 
The outer halo is fit with a single absorbed APEO model. 

boxes. This is because zo is highly sensitive to the mean 
surface brightness in the innermost bins, which is overes- 
timated according to these corrections. Overall, a scale 
height of zo > 2 arcmin (~ 4 kpc) is likely in the outer 
halo. The density scale height is t wice this value, or 
H > 8 kpc, which is consistent with Brcgm an fc Pildisl 
( 1994) considering they included all diffuse emission. In 
any case, a scale height below 1 arcmin is strongly disfa- 
vored. 

Finally, we note that there is a patch of diffuse emis- 
sion to the northwest of the galaxy with approximately 
a 12 arcmin separation. This is visible in both the Chan- 
dra and XMM-Newton data and the emission is centered 
near a =02h21m32.2s, S = +42d25ml2s. The emis- 
sion is roughly circular with a radius of about 2 arcmin, 
but it has some substructure with a brighter "bar" run- 
ning roughly parallel to NGC891. The emission is near 
the Chandra point source CXOSEXSIJ022142.6+422654 
(which may also be 1AXG J022139+4226, a source iden- 
tified as 0.5 arcmin away but in an error circle of 
1 . 1 arcmin) . We have not found any records in the liter- 
ature of this source despite its appearance in the ACIS-I 
chips of the archival Chandra data, and unfortunately 
the optical and IR data at this position is limited. In 
2MASS and Palomar/DSS images, we see no obvious 
signs of clustering. An Xspec APEC fit with only Galactic 
absorption finds a good fit for kT ~ 3keV, so we spec- 
ulate that this is a galaxy cluster at unknown redshift. 
The 0.3- lOkeV flux is F x ~ 1.3 x 10" 13 ergscm- 2 . 

2.4. Source Spectra 

We extracted inner and outer halo spectra from 
the regions described above in point-source subtracted, 
exposure-corrected data. The Chandra spectra were sep- 
arately extracted from obsIDs 794 and 4613 and com- 
bined, with corresponding response files generated. We 
likewise extracted background spectra from regions near 
the edge of the ACIS-S3 chip. The XMM-Newton spectra 
were extracted from analogous regions, and we generated 
their response files using the XMM-ESAS scripts mos- 
spectra, pn-spectra, mos_back, and pnJjack. All of the 



spectra were binned to 25 counts per bin for fitting. We 
restrict the fitting bandpass in the Chandra data from 
0.4— 2.0 keV. Above 2.0 keV, there is little source flux and 
the data do not constrain any fit. Likewise, the XMM- 
Newton data are fitted between 0.3 — 3.0 keV in the MOS 
and 0.4 — 3.0 keV in the pn. The lower bounds are deter- 
mined by the rise of the low-energy particle background 
in each detector. 

We use Xspec vl2.7.1 for spectral fitting ()Arnaudl 
[19961) . We use x 2 as the goodness-of-fit parameter, and 
define an acceptable model as one that is not ruled out 
at 95% confidence based on the x 2 value. Errors on the 
model parameters are determined using the Xspec step- 
par task, and model fluxes are quoted for a 0.3 — 3keV 
bandpass. 

The diffuse source emission is made up of two or more 
components: the bulge and (inner/outer) halo. We be- 
gin with the simplest case of a homogenous halo, but 
the inner and outer halo may be produced by a galac- 
tic fountain and IGM accretion respectively, and so may 
have different properties. The halo has an X-ray spec- 
trum that peaks near 0.6 keV, whereas the hotter bulge 
has a spectrum peaking near 0.9 keV (Figure [5]). The 
bulge emission is unimportant in the outer halo, but con- 
tributes near the disk. As we expect these to be thermal 
plasmas in collisional io nization equilibriu m, our default 
model is APEC/vapec (jSmith et alJl2001[ ). Since photo- 
electric absorption is important, as an abso rption model 
we ad opt TBabs with the abundances of iWilms et al.l 
(2000). All fits incorporate the Galactic absorption col- 
umn of N H = 6.5 x 10 20 cm- 2 (jKalberla et all [20051 ) in 
addition to variable absorption near the disk. The influ- 
ence of these choices (by comparison with alternatives) 
is described in Section [3l 

The ULX in the XMM-Newton data contributes ap- 
proximately 55,000 counts in the pn and 17,000 in each 
MOS camera in the 0.3 — lOkeV bandpass, peaking near 
1 keV. Although we mask the ULX to the 90% encir- 
cled energy radius, even outside this it contributes thou- 
sands of counts. We ad opt as a ULX model t he be st fit- 
ting simple model from IHodges-Kluck et al.l (|2012D : ab- 
sorbed emission from a hot accretion disk. This model, 
TBabs*DISKPBB, is an excellent fit to the ULX spec- 
trum (x 2 = 1169.1/1170 d.o.f. for the combined MOS 
and pn spectra) . We freeze the model parameters to their 
best-fit values: N H = 2.3 x 10 21 cm" 2 , T in = 1.62 keV, 
and p = 0.54 (p is a free parameter and T(r) oc r~ p ; 
iMineshige et al.l 119941 ) . This model is included in all 
our XMM-Newton fits with the normalization allowed to 
vary, but we find that it reproduces the expected con- 
tribution of the ULX flux very well in annuli around 
the ULX. The ULX spectrum does vary slightly with ra- 
dius because the PSF is energy dependent. However, for 
the relevant energies, the expected difference is approxi- 
mately 1%, which is smaller than the Poisson uncertainty 
in the ULX flux in the source extraction regions. 

2.5. X-ray Background 

The X-ray background includes instrumental noise, 
particle backgrounds, the cosmic X-ray background, the 
contribution from the Local Bubble, and the hot Galactic 
halo. The latter two are particularly important for our 
purposes because the Galactic halo has a similar temper- 
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Fig. 6. — XMM-Newton and Chandra spectra of the outer halo (Left) and inner halo (Right), with the MOS 1+2 in black and red, pn 
in green, and Chandra ACIS in blue. The best fits are overlaid with their residuals below, and the data have been rebinned to 5<r for 
plotting purposes. Note that the Chandra inner halo spectrum contains much more bulge emission than the XMM-Newton spectra where 
it is mostly excised with the brightest ULX emission. 



ature (kT ~ 0.25 keV) to that of NGC891, and the Lo- 
cal Bubble contributes to the emission lines by which we 
determine the temperature of the source plasma. Such 
emission is typically subtracted by using a scaled spec- 
trum extracted from nearby source-free regions. How- 
ever, the emission from NGC891's halo covers most of 
the ACIS-S3 chip and the central chips of the MOS 
detectors, and the true scale height of the halo is un- 
known. The XMM-Newton particle background also in- 
duces strong fluorescent lines (the Al and Si Ka lines 
at 1.49 keV and 1.75 keV are the important lines in our 
fitting bandpass) that are highly variable with position 
such that scaled background subtraction leaves large 
P Cygni type profiles in the 1 — 2 keV range. Thus, the 
background regions for the Chandra data contain an un- 
known amount of source flux, while those for the MOS 
detectors (placed farther out) do not directly subtract 
the instrumental background for the chip of interest. 
Our analysis fits use the background-subtracted spectra, 
but we perform parallel fits using background subtracted 
spectra and raw spectra where most of the background 
is fit explicitly (with the unexposed chip backgrounds 
subtracted, follo wing for the Chandra and XMM-Newton 
data respectively iMarkevitch et al. 2003; Snowdc n et al.l 
l200l . 

In these fits, we explicitly model the astrophysical 
backgrounds following the XMM-ESAS recipes. We use 



the HEASARC X-ray Background Toofl to estimate the 
background surface brightness from the ROSAT All- Sky 
Survey (RASS) data. This tool also provides a spectrum 
and response file, which we use to obtain the flux in each 
component of the background. 

With this tool, one usually defines a large region 
around the object of interest, assuming that the back- 
ground is approximately correct for the line of sight. 
However, there are a few bright diffuse sources within 
30arcmin of NGC891 (most notably the galaxy cluster 
containing 3C66B to the north). Thus, we measure the 
background in a circle with radius 30 arcmin from an 
off-position to the southwest. We use the flux in the 
R45 and R67 bands (0.44 - 2.04keV). We convert the 
surface brightness to energy units via the HEASARC 
WebPIMMS too0 with a 0.1 keV APEC model for the Lo- 
cal Bubble, a 0.25 keV absorbed APEC model for the hot 
Galactic halo, and an absorbed power law (r = 1.46) 
for the cosmic X-ray background. This yields, in the 
0.3 — 3.0keV bandpass, a surface brightness oil = 4±lx 
10 -15 ergs -1 cm -2 arcmin -2 . As a consistency check, we 
measured the surface brightness of the background in a 
targeted 31.5 ks ROSAT PSPC o bserva tion of NGC891 
published by iBregman fc Pildisl (|1994f ). We measured 
the 0.1 — 2.4keV count rate in a total of 2507T arcmin -2 in 
the exposure-corrected field, avoiding point sources and 

4 http:/ /heasarc. gsfc.nasa.gov/docs/heasarc/xrayback. html 

5 http:/ /heasarc. gsfc.nasa.gov/TOols/w3pimms_pro. html 
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TABLE 3 
1-T Halo Model Fits 



Spectrum 


N H ,i 
(10 20 cm" 2 ) 


fcTi Z\ 
(kcV) (Z ) 


N H ,2 

(10 20 cm" 2 ) 


kT 2 
(keV) 


z 2 

(Zq) 




X 






Outer Halo ( 


z > 2.8 kpc) 








Chandra 
Airnc i _i_o 

pn 
Joint 


6.5(f) 
b.b(t) 
6.5(f) 
6.5(f) 


n 24+ 07 04+ - 05 

u,z ^-0.03 u ' u4 -0.02 

n ni ± n no fl 1 1 +0.05 
U.zl it U.UZ U.ll q 

no-, +0.03 n na +0.05 

U - 2i -0.02 U ' US -0.03 

0.21 ±0.02 0.08 ±0.03 








102.6 (105) 

104. (^loUJ 

166.1 (175) 

426.2 (460) 






Inner. Halo (0.6 


< z < 2.8 kpc) 








Chandra 
MOS 1+2 
pn 
Joint 


10+f 
24 +27 
z ^-18 
io+22 
iS -12 

lel^ 6 


0.22 ±0.02 0.12lg;Jg 

0.19 ±0.05 o.io±g;g 7 

0.19 ±0.04 0.11±p o? 
0.21 ±0.02 O.lll^Jg 


77+12 
' '-10 

77(f) 

77(f) 
70+11 


n ca+0.06 

U - b8 -0.03 

0.58 (f) 
0.58 (f) 
0.58 ±0.05 


1.0(f) 
1.0(f) 
1.0(f) 
1.0(f) 


74.0 (79) 
155.5 (152) 
117.2 (97) 
342.15 (328) 



Note. — The MOS and pn fits include a TBabs*diskpbb component with all parameters except flux frozen as 
described in the text. The fits corresponding to these models arc shown in Figure [6] Errors are quoted at the 
90% confidence interval based on the Xspec task steppar, and (f ) designates a frozen parameter. 



bright emission to the north and west. This yields a sur- 
face brightness J = 3 ± 1 X 10 -15 ergs -1 cm -2 arcmin -2 
in the 0.3 — 3keV band, which is in agreement with the 
RASS value. 

With the normalization in hand, we fitted the RASS 
spectrum with three components: the Local Bubble (a 
thermal component with kT = O.fkcV), the Galactic 
halo (an absorbed APEC component with kT = 0.25 keV), 
and the cosmic X-ray background (an absorbed power 
law with r = f.46). Although these values may not 
be exactly correct, there are only a few spectral bins, 
so we freeze them at the values quoted in the XMM- 
ESAS manual. We use the Galactic absorbing column of 
Nh ~ 6.5 x 10 20 cm -2 for the Galactic halo component. 
In the 0.3 — 3keV bandpass, the fit finds that the Lo- 
cal Bubble contributes ~6% of the flux, the (absorbed) 
Galactic halo ~34%, and the cosmic X-ray background 
~56%, which we use to determine the normalizations of 
the model components in the source regions. 

To verify that the Chandra and XMM-Newton back- 
ground spectra do not have unknown components, we 
fitted the same model (adding Gaussian components for 
Al and Si Ka) but allowed the parameters to float. We 
find generally good agreement, but the fractional con- 
tribution of the Local Bubble is poorly constrained be- 
cause the spectrum of the kT — 0.1 keV plasma drops 
rapidly above 0.3 keV. However, the 0.3 — 3keV sur- 
face brightnesses are in agreement with the ROSAT val- 
ues: S « 4.4 x 10 _15 ergs -1 cm' 2 arcmin -2 in the MOS, 
S ~ 3.7 x 10 -15 ergs -1 cm -2 arcmin -2 in the pn, and 
S ~ 3.9x 10 -15 ergs -1 cm -2 arcmin -2 in Chandra, where 
these model fluxes are quoted for 0.3 — 3keV but fitted in 
the bandpasses noted above. These results suggest that 
the NGC891 halo surface brightness does indeed drop 
below the background on the scale of the ACIS-S3 chip 
and that the regions used for background subtraction 
are appropriate. Any overestimate of the background 
flux must be slight, and there do not appear to be large 
discrepancies between the detectors. 

Since the RASS background surface brightness and 
spectrum appear consistent with the other data, for our 
explicit background fits we fix the total flux, scaled to the 



source region, in the X-ray background. In Section [31 we 
show that fits with background subtraction and fits ex- 
plicitly fitting the background are in good agreement. 

3. SPECTRAL FITS & HALO METALLICITY 

In this section we present our model fits. Unless stated 
otherwise, we refer to the background-subtracted spec- 
tra. First, we describe single-temperature (1-T) fits to 
the outer and inner halo, followed by a more complex 
two-temperature (2-T) model. We then compare the re- 
sults of our 1-T fits to the fits obtained by fitting the 
background explicitly. Finally, we discuss several po- 
tential sources of systematic error and the effect on our 
results. 

3.1. Models 

The models we use for the different spectral compo- 
nents are described in the previous section and summa- 
rized here. 

The outer halo is represented by an absorbed ther- 
mal (TBabs*apec) model, with N H fixed at the Galac- 
tic value. This model represents emission from a hot 
plasma in collisional ionization equilibrium modified by 
photoelectric absorption by an intervening column of 
atomic gas. The inner halo is the same, but with Nh 
variable. The bulge uses the same model, but with 
Z = Zq. For the ULX, we use a n absorbed accretion 
disk m odel (TBabs*DISKPBB) as in lHodges-Kluck et"aL 
( 2012f) and described above. Following iSnowden et al. 
(2004), we model the Local Bubble with an APEC model 
of kT = 0.1 keV. The Galactic halo and cosmic X-ray 
background are modeled as TBabs(apec+pow), where 
Nh is set to the Galactic value, kT — 0.25 keV, and 
T = 1.46. The total flux is scaled to the RASS value, 
and the power in each component is constrained by sim- 
ulaneously fitting the RASS spectrum. The Al and Si 
Ka lines are fit by GAUSS components, with fixed energy 
and a — (instrumental linewidth). 

Thus, the 1-T model to the Chandra outer halo 
is just TBabs*APEC, whereas in XMM-Newton data 
it is TBabs*APEC+TBabs*DISKPBB. In the inner 
halo, we add the bulge component. In the fits 
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Fig. 7. — Left: Z vs. kT contours for the outer halo (top) and inner halo (bottom). The contours represent 1, 2, and 3cr deviations from 
the best-fit values denoted by crosses. The shaded regions show the contours for the joint fit. Right: As left, for the flux (normalization) 
vs. Z contours. The normalizations have been scaled to the Chandra value (due to different areas) and the absolute value is unimportant. 



where the X-ray background is explicitly fitted, we add 
APEC+TBABS(APEC-r-POW) components and Gaussians 
in the XMM-Newton data. The best fit values for the 
background-subtracted fits are tabulated in Table [3] 

3.2. Outer Halo 
The Chandra spectrum (Figure [6]) is fit well by a 
thermal model with kT = 0.24t°;^keV and Z = 

°- O4 ^ao2^0- The M0S and P n spectra prefer a slightly 
cooler (kT ~ 0.21 keV; Table [3]) plasma with a metallic- 
ity Z ~ O.IZq, but the results agree within their 90% 
error bars (for a single parameter of interest). This is 
cooler than the . 3 keV found in the ROSAT data by 
iBregman fe Pildii ((1994T ) . although most of the light in 
their spectrum comes from the inner halo region. Like- 
wise, it is slightly c ooler than the 0- 26 keV XMM-Newton 
value reported by iTemple et al.l ()2005l ) using a much 
shorter exposure, but again their spectrum is dominated 
by the inner halo. All of the fits are excellent (Figure [6]) 
and the temperatures are constrained to 15%. 

The higher signal and higher spectral resolution of the 
XMM-Newton data also mitigates the degeneracy be- 
tween kT and Z in the Chandra data, as seen in a 2- 
parameter A% 2 contour plot (Figure [7]) . It is evident 
from this plot that solar metallicity is ruled out at more 
than 3cr confidence in each detector. There is an impor- 
tant degeneracy between the flux and metallicity (the 
normalization values have been scaled to the Chandra 
value due to different extraction areas, and the abso- 



lute values are unimportant), but subsolar metallicity is 
preferred as strongly as in the Z — kT plot, where the 
degeneracy between the parameters (which manifests as 
non-zero slope) is much smaller. Since the results from 
each detector are largely in agreement (Table |3]) , we can 
jointly fit all four spectra with a 1-T model. With a total 
of 24602 source counts in the fitting bandpass (of which 
we estimate fewer than ^6000 arc ULX photons), we find 
5a limits of kT = 0.2Q±g;g| keV and Z = 0.11±J];^Z Q for 
a single parameter of interest. 

We measure the relative abundances of O and Fe by us- 
ing a thermal model where the abundances of each metal 
are variable (vapec). Since all of the elements contribute 
to the continuum while only some have strong lines in the 
soft X-ray bandpass, we allow only O and Fe to vary, ty- 
ing the abundances of like elements to their values and 
setting the rest to the best- fit Z value in Table [3] for each 
spectrum. We note that if the other elemental abun- 
dances are tied together but allowed to vary, the best-fit 
O/Fe value is very similar to that obtained when these 
abundances are frozen. 

The MOS and pn spectra indicate O/Fe~0.7, whereas 
the Chandra spectrum is consistent with Ap c = 0, so 
O/Fe is greater than unity but poorly constrained. The 
fit contours are shown in Figure |SJ As in the APEC fits 
(Figurc[7]), the MOS has a higher overall metallicity than 
the pn, but the ratio of O/Fe is similar between the two 
(0.69 in the MOS and 0.67 in the pn, well within the 
statistical error) . The Chandra spectrum prefers a lower 
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TABLE 4 
2-T Model Fits 



Spectrum 


fcTi 
(keV) 


kT 2 Absorbed 0.3-3.0 keV 
(keV) F 2 /F 1 


Norm. Ratio 
EM 2 /EMi 


x 2 


Outer Halo (z > 2.8 kpc) 


Chandra 
MOS 1+2 
pn 
Joint 


09+ 06 
u,ua -0.02 

12+ ' 04 

0.11 ± 0.02 

0.11 ± 0.02 


0.29 ±0.05 1.9 ±0.7 
0.25±g"g| 3.4 ±0.3 
0.27 ± 0.03 2.5 ±0.3 
0.27 ±0.03 - 


0.33 ±0.07 
0.29 ±0.06 
0.37 ±0.07 


109.4 (105) 
173.2 (180) 
169.1 (175) 

436.5 (460) 



Inner. Halo (0.6 < z < 2.8 kpc) 



Chandra 0.09^°;^ 



MOS 1- 
pn 

Joint 



2 0.13 ±0.05 
■0.06 
0.03 
,+0.05 
0.03 



0.05] 
0.09j 



0.25 ±0.05 5±2 

0.26 ±0.06 2.6 ±0.4 

0.21 ±0.05 3.1 ±0.7 

0.24 ± 0.03 - 



0.3 ±0.1 
0.5 ± 0.1 
0.4 ±0.2 



72.2 (79) 

152.3 (152) 
111.9 (97) 

346.4 (328) 



Note. — The fits include components as in Table[3] including the ULX in the XMM-Newton data 
and a thermal bulge component in the inner halo. Errors are quoted at the 90% confidence interval 
for one parameter of interest based on the Xspec task steppar, and for each kT value the other is 
frozen during the error search to prevent reversals. Column 4 is the 0.3 — 3 keV absorbed flux ratio 
between the two thermal components, and column 5 is the ratio between their model normalizations. 
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Fig. 8. — 2-parameter contour plots for O/Fe in the outer halo 
(Left) and inner halo (Right). The contours are 1, 2, and 3cr con- 
tours and the minima are marked by crosses. The shaded region is 
the 1, 2, and 3cr contours for the joint fits. 



Fe abundance because of the rather low flux between 
0.7 — 1.0 keV (Figure [6]), but the value is poorly con- 
strained, so the joint fit settles on a value near O/Fe~0.7 
(Figure [5]). However, O/Fe is quite sensitive to tem- 
perature (higher temperatures tend to depress Ap c and 
increase Ao ) . The error bars should therefore be consid- 
ered to be somewhat larger than the contours in Figured! 

3.3. Inner Halo 

Despite its higher surface brightness, the inner halo 
contributes far fewer source counts due to the smaller 
regions we have demarcated. The XMM-Newton spec- 
trum is also more contaminated by the ULX and both 
the Chandra and XMM-Newton spectra include contam- 
ination from the bulge emission (seen especially in the 
Chandra spectrum in Figure [5]). Contamination from 
the disk is minimal because of our extraction region 
boundaries (outside two disk scale heights) and the very 



high absorption towa rds the disk (Nh ~ 10 cm in 
lOosterloo et~aill2007l) . Finally, N H varies across the re- 
gion, so it must remain a free parameter in spectral fits. 

We adopt a model analogous to the outer halo where 
we include a bulge component and allow the Nh column 
to vary above the Galactic value. As the bulge is nearly 
cospatial with the ULX in the XMM-Newton data, most 
of its emission is masked along with the ULX. The re- 
maining bulge contribution is dominated by the Fe L- 
shcll lines between 0.7 — l.OkeV, which can be fit by a 
range of models with high absorption and moderate tem- 
perature. We therefore constrain the bulge temperature 
using the Chandra data. A 1-T model is a poor fit to the 
Chandra data, so we use a 2-T fit with the bulge metallic- 
ity frozen at solar (the metallicity is poorly constrained, 
but consistent with solar and well above the outer halo 
value). This yields kT ~ 0.6 keV (TableEJ). For the bulge 
contribution in the MOS and pn, we use a thermal com- 
ponent with the absorption and temperature frozen at 
the Chandra values and the flux free to vary. However, 
fits with unrealistic bulge fluxes (i.e., comparable to or 
more than the total Chandra bulge flux) are rejected. 

Using this model, we find good fits for thermal mod- 
els in the EPIC spectra that are consistent with the 
Chandra data (Table [3]). These fits and their residu- 
als are shown in Figure [6j The best- fit temperatures are 
kT m 0.21 keV, while the metallicity is near Z ~ 0AZ@ 
in each spectrum. Because of the lower signal, the con- 
fidence intervals for these parameters are substantially 
wider than in the outer halo (Figure [7]). We note that it 
is possible to find good fits with different combinations 
of kT, Z, and Njj with only small differences in % 2 , so 
while the 2-parameter contours are accurate for the best 
fit, the 3-parameter region of parameter space with good 
fits is substantially wider. As might be expected, the 
temperature is known much better than the metallicity. 
Still, it is noteworthy that the best-fit parameters in each 
detector agree not only with each other (within the 90% 
error bars) but also with the outer halo values (Table |3]). 

As in the outer halo, we can improve the signal by 
jointly fitting the four spectra. This yields 3tr limits of 
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Fig. 9.— Top: Folded EPIC MOS 1+2 outer halo spectrum fit 
by a 2-T galactic fountain model with components as tabulated 
in Table [4] The dashed line shows the contribution of the ULX. 
Bottom: As above, for the 1-T model shown in Table [3] 



kT = 0.21 



+0.05 



keV and Z = 0.11" 



The absorb- 



-0.06 ^ v ^ — "• ii -o.o7 Zy 0- 
ing column is poorly constrained even in the joint fit, 

with a 3cr limit of N H = 16+ 26 x 10 20 cm" 2 on top of 
the Galactic value. Based on the H I data (cf. Figure 
we expect a value above Nh = 10 21 cm' 2 but below 
10 22 cm -2 , although generally we expect more emission 
to come from less absorbed regions, biasing the fitted Nh 
value toward lower values. Thus, we cannot rule out a 
higher metallicity for the inner halo than the outer halo, 
but both are likely to be subsolar. 

The signal is not sufficient to reliably measure the rel- 
ative abundances if the bulge component is also allowed 
to vary. By freezing the bulge component at the best- 
fit apec values, we can obtain O/Fe values, but these 
should be interpreted with caution. We proceed as with 
the outer halo, and Figure [8] shows the resultant con- 
tours. Like in the outer halo, the best fits tend to have 
0/Fe< 1.0, although the pn now finds a value of O/Fe 
slightly smaller than 1.0. The Chandra spectrum indi- 
cates a higher Fe abundance than in the outer halo, in 
agreement with the APEC fit (Table |3]). However, the 
uncertainty is large and the bulge emission is quite im- 
portant to this spectrum. If we fit the total Chandra 
spectrum extracted from both the inner and outer halo, 
the best-fit value is then 0/Fe= 1.3lg 4- The la contour 
for the joint fit is in agreement with the joint fit for the 



Fig. 10. — Simulated 400 ks XMM-Newton spectra created from 
the best-fitting 2-T model (excluding the ULX). The simulated 
spectrum is fitted with the 2-T (solid line) and 1-T (dashed line 
and dotted line) models, with the residuals for each below. The 
middle residual panel shows a 1-T fit with Z frozen at 0.12Zq (as 
in Table |3j , and the bottom panel has variable Z with a best-fit 
value of Z = 0.28.Zq. Note that the i/-axis is linear rather than 
logarithmic, and that the data have been binned to 10rj significance 
for plotting purposes. 



outer halo, but the best fit finds a slightly higher O/Fe. 
The picture is confused by the bulge emission, but at face 
value the O/Fe value is slightly higher in the inner halo 
than the outer halo. 

The best-fit values suggest that the inner and outer 
halo are composed of the same plasma despite the ap- 
parent connection between the inner halo morphology 
and the activity in the disk. "Deprojecting" the Chan- 
dra spectrum by using a scaled outer halo spectrum as 
the background confirms that the inner halo tempera- 
ture we derive is not merely due to outer halo material 
in the line of sight. Thus, the hot halo of NGC 891 is con- 
sistent with being an isothermal, 2 — 3 x 10 6 K plasma. 
The agreement with the outer halo indicat e s tha t the 
temperatures report ed by iBregman fe Pildisl (|1994l) and 
iTemple et all (|2005D of 0.3 keV and 0.26 keV respectively 
may be boosted by contamination from the bulge. How- 
ever, they are consistent with the limits obtained by fit- 
ting the spectra individually, and there is some degener- 
acy between Nh, kT, and Z. 

3.4. A 2-T Model 
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Fig. 11.— Top: Folded EPIC MOS 1+2 outer halo spectrum 
with background components explicitly fitted and residuals. The 
dotted black line is the APEC model for NGC891's halo, while the 
dashed black line is the ULX contribution. The blue lines are 
background components: the Local Bubble (dash and dots), the 
Galactic halo (dashed), the cosmic X-ray background (solid), and 
the instrumental fluorescent lines (dash-dot). The power in the X- 
ray background components is constrained by the fit to the RASS 
spectrum. Bottom: As above, but for the background-subtracted 
(analysis) spectrum. 



It is also possible to obtain a good fit with a 2- 
T model where the metallicity of both components is 
fixed at solar (the bulge emission is treated as above). 
This fit is motivated by the galactic fountain scenario as 
well as the many 2-T fits in the liter ature to the halos 
(|Strickland et al.ll200l and disks (e.g. iTvler et al.ll200l 
of massive late-type galaxies. The fit requires the first 
component to have kT\ ~ 0.1 keV and the second to have 
kT2 ~ 0.25 keV (parameters given in Table |4j. The ratio 
of the emission measure is weighted towards the cooler 
component. 

As in the 1-T model, the inner and outer halo spectra 
have parameters consistent with each other. The fits 
in Table 0] are carried out in the same way as in the 1-T 
model except that we add a second APEC component and 
fix the ULX flux in the XMM-Newton spectra to the best- 
fit value in the 1-T model. The error on each temperature 
is computed while freezing the other because they are 
separated by only 0.15keV and can exchange places in 
the parameter space search. 

As the 1-T and 2-T models both have well constrained 
parameters and produce good fits, it is worth asking how 
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Fig. 12. — Comparison of Z — kT 1, 2, and 3<r contours for the 
background-subtracted and background-fitted outer halo spectra 
for the MOS (Left) and pn (Right). The solid black line represents 
contours from the background-fitted spectra and the dashed red 
line the background-subtracted (analysis) spectra. Crosses identify 
the minima. 



much complexity the spectrum can distinguish and how 
much signal is required to rule out one of these spectro- 
scopically. To examine this question, we used the Xspec 
tool fakeit to simulate spectra based on the best-fit 2-T 
model in the outer halo using the pn response, increas- 
ing the signal by increasing the exposure time. In these 
simulations, we did not include the ULX. We created 10 
such spectra for each exposure time, then fit them using 
the 1-T model. 

We find that twice the S/N (-400 ks) is enough to 
distinguish between the 1-T and 2-T models to high con- 
fidence. The difference between the models can be seen 
in the residuals (Figure [TO)) , where the 1-T model un- 
derpredicts the O VII flux and overpredicts the O VIII 
flux. Above 1.0 keV, the models have similar residuals, 
and if the ULX were included it would dominate the 
spectrum here. Even at lower signal, it may be possi- 
ble to rule in favor of one model based on external cri- 
teria or observable predictions of each model. We ex- 
amine a few of these in Se ction |4j We also note, as in 
IStrickland fc Stevens! (|2000( ) , that the halo spectrum may 
not be simple and could be undermodeled in either case. 

The existing data are good enough to rule out non- 
thermal models (i.e., those that do not produce emission 
lines). The structure in the spectrum below 1.0 keV can- 
not be reproduced by realistic absorption of a smooth 
continuum such as a power law. Continuum models pro- 
duce poor fits with residuals similar to, but worse than, 
those seen in Figure LTOl 

3.5. Background Fits, X-ray Binaries, and Possible 
O VII Excess 

For the fits described above, we have used the 
background-subtracted spectra. To verify these fits, we 
also fitted the raw spectra (with model quiescent par- 
ticle backgrounds subtracted) while explicitly fitting the 
X-ray background, follo wing the XMM-ESAS guide (and 
iMarkevitch et a l. 2003, for the Chandra data) and using 
the RASS data. We also check the background fit to the 
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Fig. 13. — Left: 1, 2, and 3er Z — kT contours for fits to the 
XMM-Newton and Chandra outer halo spectra as reported in Ta- 
ble [3] Right: As left, but the XMM-Newton fits incorporate a 
Gaussian component at 0.58 keV representing the O VII complex. 
The Chandra contours are the same as at left. See text for discus- 
sion. 



source regions against background spectra extracted for 
each detector. 

The results in the outer halo for the MOS are shown 
in Figure 111! In the top panel, the background compo- 
nents (folded through the MOS response) are overplotted 
in blue while the source components are shown in black. 
The spectrum is fitted simultaneously with the RASS 
spectrum with the non-flux parameters tied together to 
constrain the power in each component. The overall 
flux in the background is computed from the surface 
brightness measurements in the RASS and MOS back- 
ground, so if the RASS model (Local Bubble, Galactic 
halo, and cosmic X-ray background) is a good description 
of the background, the source fits should be identical to 
the background-subtracted case. FigurelTTIdemonstrates 
that both cases give excellent fits, and Figure [12] shows 
that both cases essentially give the same values for Z 
and kT. This is further verified by fitting the MOS and 
pn background spectra, which are fit well by the RASS 
model. 

The effect of the Local Bubble or Galactic halo on 
the measured temperatures and metallicities must there- 
fore be small, as must be the contribution of the halo 
of NGC 891 beyond the ACIS-S3 chip or central MOS 
chips. This is particularly important considering that 
the Galactic halo has a temperature (kT ~ 0.25 keV) 
similar to NGC 891 's halo. 

We have thus far not considered the flux contributed 
by X-ray binaries (XRBs). This flux may be from un- 
detected sources or from the wings of the PSF of de- 
tected ones. Ignoring the XRB contribution can ar- 
tificially reduce the metallicity by forcing the thermal 
model to boost continuum flux (for a discussion, see 
IKim & Fabblano| [2063f). The tell-tale sign of a missing 
XRB component in a thermal spectrum is positive resid- 
uals at high energies. The residuals to our background- 
subtracted halo fits (Figure [5]) do not indicate that an 
additional component is required, except perhaps in the 
inner halo spectrum of the pn. The MOS and pn outer 



Fig. 14. — Red: 0.5-0.6 keV adaptively smoothed XMM-Newton 
image showing emission near the O VII complex. Green: 2MASS 
H-band image. Blue: Total H I column density, with contours 
overlaid at 10 19 , 10 20 , 10 21 , and 10 22 cm" 2 . 



halo spectra show a similar structure in the residuals 
above 2 keV, but this appears to be a consequence of 
the background subtraction (Figure ITT]) • As the XMM- 
Newton spectra contain a ULX contribution, the Chan- 
dra data are a better indication of the influence of XRBs. 
Here, too, there is no obvious need for a harder compo- 
nent. This is essentially what we would expect in the 
halo, since the XRB population should predominantly 
occupy the disk. 

From a model perspective, no additional component is 
required (as with the 2-T spectrum). However, to see 
if an XRB component materially improved the fit, we 
extracted an aggregate spectrum from the brightest non- 
ULX disk sources and fit it with an absorbed power-law 
model. This yields a photon index of T ~ 1.8, similar to a 
power-law fit to the ULX. We then fitted the background- 
subtracted halo spectra with a thermal model and a 
power-law component with a common absorbing column 
and r = 1.8. In the XMM-Newton spectra, the compo- 
nent is degenerate with the ULX, and in the Chandra 
spectrum the improvement to the fit is small and the 
change to Z is within the 9 0% error bars reported i n Ta- 
ble ©. Thus, in contrast to IKim fc Fabbianol (|2003l ). our 
finding of low metallicity in the halo is not likely influ- 
enced by the XRB population. Fits with T = 1.4 (similar 
to the CXB) demonstrate that the AGN contribution is 
also small. 

As mentioned above, the residuals in the pn and MOS 
fits (Figured]) suggest that the O VII flux near 0.58 keV 
is underpredicted while the O VIII flux is overpredicted. 
Indeed, adding a Gaussian component of zero (instru- 
mental) width to the XMM-Newton spectral fits does 
improve the fit (e.g. in the pn by A% 2 ~ 8 when remov- 
ing a single de gree of freedom for the flux of the line). 
As discussed in iProtassov et all (]2002l ). the F-test is in- 
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appropriate for determining the presence of a line, but 
when the parameters are well constrained (as they are in 
our model), the F-test can be "calibrated" using simu- 
lated spectra. Using the fakeit tool in Xspec, we simulate 
100 spectra based on the best-fit (null) model and pn re- 
sponse and exposure time to determine the fraction of 
data sets in which the O VII line improves the fit by a 
similar amount. There are several cases in which this 
occurs, but none have a A% 2 larger than 8. Finally, the 
O VII complex is detected in the RGS data. 

The presence of the "line" is unsurprising if the spec- 
trum does indeed require a 2-T fit. However, it is unclear 
that the emission is attributable to NGC 891. This is be- 
cause it does not appear to be necessary in the Chandra 
spectra (although these have worse spectral resolution) 
and adding the line as a Gaussian component brings the 
XMM-Newton 1-T model parameters into better agree- 
ment with the Chandra values (Figure [T3|) . In the right- 
hand plot, the line has been added and the temperatures 
agree. 

On the other hand, from a XMM-Newton combined im- 
age made in the 0.5 — 0.6 keV band (Figure [T4|) . there is 
some weak but suggestive evidence that the strong O VII 
emission is related to the galaxy. The emission basically 
tracks the inner halo (as one expects, considering the halo 
continuum is bright in this band), but there is a small ex- 
tension to the northeast following the outermost H I con- 
tour. The hole in the middle is the ULX mask. Although 
the image shown is adaptively smoothed, this statement 
is based on adaptive binning. The Chandra data do not 
show many point sources missed by the XMM-Newton 
mask in this region, so the excess is likely diffuse or made 
up of faint point sources. 

If the emission is not connected to NGC 891, it is 
also not likely from the Local Bubble or Galactic halo. 
Although the Local Bubble should radiate strongly in 
O VII, to explain the excess it would need a flux similar to 
that in the 2-T model, or 3— 4 times that predicted by the 
RASS data. Since the Galactic halo has a temperature 
similar to the 1-T fit temperature (kT ~ 0.2 - 0.25 keV), 
it likewise cannot produce a strong O VII line without 
also producing a strong O VIII line, which is the same 
reason the 1-T model leaves this excess (Figure [TDJ . 

It is probably not produced by solar wind charge 
exchange because, while we would expect charge ex- 
change to produce a simple line, we expect charge ex- 
change to produce other strong lines (e.g. O VIII), and 
ICarter fc Sembavi (|2008D find that XMM-Newton obser- 
vations affected by charge exchange have line light curves 
that are uncorrelated with continuum curves. We find 
nearly perfect correlation between the O VII line and 
nearby continuum (a Pearson coefficient of R — 0.997). 
There is also no record of instrumental O VII lines in the 
XMM-ESAS manual or calibration reports. 

Overall, it is unclear whether the O VII flux supports or 
enables the 2-T model because it seems to be peculiar to 
the XMM-Newton data and brings the 1-T temperatures 
into agreement with the Chandra data. 

3.6. Systematic Error in kT and Z 

As the metallicity is the crucial discriminant between 
the accretion and galactic fountain scenarios, it is worth 
discussing the best fit in a phenomenological sense to 
gauge the uncertainty in Z. 
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Fig. 15.— 0.3 - 3keV MOS 1+2 outer halo spectra with three 
fits overlaid. The best fit (Table [3} has kT = 0.24 keV and Z = 
0.12Zq. The dashed line represents the best fit with kT frozen 
at 0.24 keV and the metallicity frozen at Zq, whereas the dotted 
line is the best fit when kT is variable (kT = 0.20 keV) but the 
metallicity is frozen at Zq . The residuals of these cases are shown 
below. See text for discussion. 



In principle, the temperature and metallicity can be 
uniquely determined from the spectrum of a plasma 
in collisional ionization equilibrium based on the lines 
present and continuum flux. Once the temperature has 
been established from the relative strength of identify- 
ing lines, the metallicity can be determined by scaling 
to the correct flux. At the low spectral resolution of 
the CCDs, the temperature is determined (below a few 
keV) by the peak wavelength of the spectrum, since the 
lines produced at different temperatures are well sepa- 
rated in wavelength. The metallicity is determined in 
the same way as before, but at low resolution, there are 
fewer patches of (effectively) pure continuum. 

The spectrum of the hot halo of NGC 891 peaks at ap- 
proximately 0.6 keV (Figure[5]). The location of this peak 
tightly constrains the temperature, with uncertainties in 
the best-fit values of around 15%. The low metallicity 
results from the substantial flux and shallow slope be- 
low 0.5 keV where there are few strong lines and the flux 
must be provided by the continuum. Above 0.5 keV there 
are few regions of pure continuum, so it is primarily this 
region that drives the fit. Figure [TS] demonstrates this 
in the MOS 1 and 2 spectra. In this plot the solid line 
is the best fit as reported in Table |3l The dashed line is 
the best fit obtained when the temperature is frozen at 
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the value in Tableland the metallicity is frozen at solar. 
The primary effect of raising the metallicity is seen in the 
residuals below 0.5 keV, where the model is a poor fit. If 
the temperature is allowed to vary while the metallicity 
remains fixed at solar, we can see from the dotted line 
in Figure [TBI that the fit is improved but still poor below 
0.5 keV. This effect would be even stronger without the 
ULX contribution. Thus, our results depend critically on 
the flux in this region. 

There are several sources of potential bias, which can 
be divided into problems with the data and problems 
with the modeling. There are some apparent differences 
between the XMM-Newton and Chandra spectra (Fig- 
ure [7]), but the differences are small in absolute terms 
and fall within the 3<7 statistical error contours. Since 
the metallicities in all three are small, if the calibration 
below 0.5 keV is wrong, it is wrong in a similar way for all 
three spectra. The ULX is unlikely to be a source of bias 
given the tight constraints on its model parameters and 
small contribution below 0.5 keV (Figure fTTj) . The agree- 
ment in astrophysical background surface brightness and 
spectra between the RASS and three instruments we use 
indicates that the soft shoulder is not likely to be due 
to inappropriate background subtraction, for example a 
higher-than-expected Local Bubble contribution. While 
fits to the instrumental backgrounds without constraints 
from RASS do not constrain the Local Bubble emission 
well (it drops sharply above 0.3 keV), in a phenomenolog- 
ical sense it does not matter how we fit the background 
as long as it is appropriately "subtracted" in the explicit 
fits, and the agreement between the parameters (Fig- 
ure indicates that it is. 

On the modeling side, one might worry about our treat- 
ment of the photoelectric absorption, the relative abun- 
dance tables we use, the thermal model we use, and the 
(non-)modeling of the possible O VII excess. 

Photoelectric absorption predominantly affects the 
softer X-rays, thus potentially biasing kT and Z . In our 
fits to the outer halo, we fix Njj at the Galactic value of 
Njj ~ 6.5 x 10 20 cm -2 . The error in this value is small 
compared to the effect of ignoring the Nh column in the 
outer halo of NGC891, which we do. Fortunately, the 
magnitude of the shifts, seen best in the Chandra data, 
is small: about ±0.01 keV and ±0.01 dex in kT and Z 
respectively for the la errors on Nh- There appears to 
be no bias towards one direction. 

The WILM abundance ta ble is an updated version of the 
lAnders &; Grevessel (|1989f ) table (see their paper for de- 
tails). The difference betw een the tables amounts to (us- 
ing the lWilms et al.ll2000L as a default) a decrease in kT 
of less than — 0.005 keV but a decrease in Z of —0 . 04 dex 
in the MOS spectra. Using the lGrevesse fe Sauvall ()1998[ ) 
table causes a relative increase of ±0.02 keV to kT while 
decreasing Z by —0.02 dex . Ther e is almost no difference 
betwe en the iWilms et all ([2000) table and the lLoddersI 
( 2003) table, which was based on updates to solar pho- 
tosphcric abundances. The primary effect of choosing a 
different table is to decrease the metallicity. 

We use the APEC and vapec thermal models with the 
ATOMDB v.1.3 transition database. Although there are 
differences between the APEC and MEKAL codefl the ef- 
fect on the kT and Z parameters in our spectra is quite 

1 see http://www.atomdb.org/Issues/mekalspex.php 



TABLE 5 
1-T and 2-T Derived Quantities 



Model kT M 0.01 - 3.0 keV L x H M tot 
(keV) (Msyr- 1 ) (ergs' 1 ) (kpc) (M Q ) 

1- T 0.2 0.4 6 x 10 39 5 ± 2 3 x 10 s 

2- T 0.1 3 4 x 10 40 3 ± 1 1 x 10 s 
2-T 0.25 < 0.05 < 2 x 10 39 6 ± 2 < 10 8 

Note. — Derived quantities discussed in Section BTTI 

small, less than ±0.01 keV and ±0.02 dex respectively. 
We note that since the models share some calculations, 
this cross check does not account for errors in some of 
the underlying formulas or values. If there is bias here, 
its direction is unknown. 

Finally, the possible O VII excess, modeled as a line, 
does change the temperature (Figure [T3"|) , since it is an 
important diagnostic line. However, if the excess is real 
it may indicate an additional thermal component rather 
than a single line, so this change in temperature should 
be treated with caution. Nonetheless, the line does not 
appear to change the overall metallicity much, as its 
presence or absence does not change the need to model 
the flux below 0.5 keV. Adding a line component does 
change the O/Fe value, counterintuitively increasing the 
value because the temperature of the 1-T model is forced 
higher, increasing the emission in the Fe L-shell lines. In 
the fit as reported in Table [3j the low temperature re- 
quires a relatively higher Fe abundance to produce these 
lines. 

Overall, the additional uncertainty incurred by our 
choice of instruments, model, abundance table, and Nh 
column do not cast doubt on the low metallicity in the 
1-T fits to the halo. The dominant source of uncertainty 
in kT and Z appears to be the degeneracy between the 
two, which is exacerbated by the low spectral resolution 
of the CCDs. 

4. DISCUSSION 
4.1. Spectral Models 

In the prior section we found that both the inner and 
outer halo are well characterized by a 1-T thermal model 
with kT ~ 0.2 keV and Z ~ 0.1Z®. Solar metallicity is 
strongly disfavored, as are absorbed continuum models 
that do not produce emission lines. However, the spectra 
also admit a physically motivated 2-T model with solar 
metallicity where a cool kT ~ 0.1 keV component is dom- 
inant but a kT ~ 0.25 keV component is required. The 
metallicity of the 0.1 keV component must be close to 
solar (or the fit reverts to the 1-T model), but the metal- 
licity of the 0.25 keV component is poorly constrained. 
In our fits in Table SJ we froze it at solar. 

We now ask whether the data distinguish between 
these models based on the cooling rate, scale height, and 
total mass involved in each model. For the purposes of 
this discussion, we consider the 0.25 keV component in 
the 2-T model to be an accreted corona (possibly with 
subsolar metallicity). 

It is worth noting that neither model may be appro- 
priate if the X-ray emission is more complex. If, for 
example, the X-rays originate in cold material that is 
shocked by a much hotter medium (as in galactic winds; 
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iStrickland fc Stevendl2000h . the 1-T and 2-T models are 
insufficient to describe the true spectrum. From a his- 
torical standpoint, there is some reason to believe this 
might be true: early, low-resolution and low-signal X-ray 
spectra from galactic halos were fit by single-temperature 
models with low metallicity, but later fits incorporat- 
ing foreground and background components relaxed the 
metallicity constraints and fou nd a 1-T vs. 2-T ambi- 
guity as we find in NGC891. IStrickland et al.l (|2004f ) 
found that the X-ray spectra extracted from a sample 
of nearby, edge-on, star-forming galaxies could all be fit 
with the same 2-T model (with different normalizations). 

However, the simplest versions of the galactic foun- 
tain and IGM accretion scenarios both predict a simple 
halo. This warrants examining the viability of these sim- 
ple models in view of quantities derived from the data. 

4.1.1. Cooling Rate 

Without heating terms, a simple estimate of M is the 
luminosity divided by the specific energ y, M = L/e, 
where e = 3/cT/2/imH (Nul sen et al.1119841 and the mean 
particle mass fx = 0.63 for a primordial composition. 
When cooling is dominated by line emission, as is the 
case for kT ~ 0.1 — 0.3 ke V, the cooling rate does no t 
depend strongly on density ( Sutherla nd fc Dopital fl993). 
As both the 1-T and 2-T models radiate substantial en- 
ergy below the 0.3 keV cutoff the Chandra and XMM- 
Newton CCDs, we use unabsorbed model luminosities 
between 0.01 — 3keV, summing values from the inner 
and outer halo. This is L ~ 6 x 10 39 ergs _1 in the 1-T 
model and L ~ 4 x 10 40 ergs _1 in the 0.1 keV compo- 
nent of the 2-T model. For reference, the 0.3 — 3keV 
model luminosity of the h alo is Ly w 2 x 1 39 er gs~ 1 , 
which is in agreement with IStrickland et al.1 (|2004f ) , and 
the 0.1-2keV luminosity L x ~ 4x 10 39 ergs _1 in the 1- 
T model is in agreement wi th the ROSAT data reported 
in lBregman fc Pildisl (|1994f ). 

The 0.01—3 keV luminosities yield cooling rates of M ~ 
0.4M Q yr _1 in the 1-T model and M ~ 3M© yr" 1 in the 
0.1 keV component of the 2-T model. The 0.25 keV com- 
ponent of the 2-T model contributes M < 0.05Mq yr _1 . 
We have corroborated these values by measuring the den- 
sity in the inner halo derived from spectra extracted near 
the disk in thin strips. We obtain the density by fit- 
ting the spectrum to yield the model emission measure, 
which we use to derive the density by assuming a vol- 
ume. The volume we use is a rectangular prism with 
two axes given by the extraction box on the sky while 
the third axis is estimated by assuming that the extent 
of the halo emission along the line of sight is the same 
as the extent of the galactic disk, which we compute as 
a chord across the circle at the mean galactocentric ra- 
dius of the extraction box. Since the emission measure 
is proportional to n 2 , the error along the third axis en- 
ters as the square root of the chord length. In practice, 
the sky axes we use are 5kpc by 2kpc, whereas we de- 
rived densities at several radii (with imputed box depths 
between 3 — 6 kpc) and found consistency within the sta- 
tistical errors. Wit h the density and cooling fun ction (we 
use the function in Mathews & Bre gmanlll978f) . we can 
compute i coo i and the total mass, thereby producing an 
average cooling rate. The values are consistent with the 
luminosity-based cooling rate, with M < 0.5M Q in the 



1-T model and M ~ 2M Q in the 0.1 keV component of 
the 2-T model (the density tends to be overestimated 
based on uncertainty in the metallicity and Njj, but we 
fix Z = Z Q for the 2-T model). 

These values are upper bounds if there are heating 
terms, in which case L becomes L not . The most plausi- 
ble so urce is hot, rising fountain material. Melio li et al.l 
(2009) find that such heating suppresses cooling in the 
accreted corona entirely. However, the 0.1 keV fountain 
material is cooler than the 0.25 keV accreted compo- 
nent, besides which it is injected trans- or subsonically 
(jHouck fe Bregmanl 119901 ) . so we expect it to cool too 
quickly to heat the wider corona. In the accretion sce- 
nario, the gas is heated to the virial temperature at the 
accretion shock, but in the steady state this cannot, by 
definition, prevent cooling. Thus, we consider the heat- 
ing to be a small effect. 

Coronal gas will cool through lines such as O VI 
AA1032,1038A, and the 0.1 keV component of the 2-T 
model is expected to produce an O VI A1032A luminos- 
ity of about L1032 ~ 7 x 10 38 ergs _1 . The cooling rate 
seen in these lines may therefore favor one model. 

Unfortunat ely, there are n o O VI detections in 
NGC891, but lOtte et all (|2003D report an upper limit of 
£1032 < 7.5 x 10 38 ergs -1 assuming that the O VI halo is 
coextensive with the X-ray halo. Since they account only 
for foreground extinction, we compute more accurate 
limits by using the average H I column densities at the 
positions of their spectra. Assuming the typical Galac- 
tic value for reddening, R(V) = A{V)/E{B - V) = 3.2, 
holds in NGC 891, the maximu m total neutral column 
seen within the lOtte et~a l. (2003) positions is E(B-V) ~ 
0.4 very near the disk and less than 0.1 farther out. 
Most of the flux in their positions comes from regions 
with between 1 — 5 x 10 20 cm -2 of intrinsic absorption. 
We extrapola te the reddeni n g val ues to the UV using 
the curves of ICardelli et~aT1 (fl989h . Using the average 
intensities and a column averaged over this region, we 
find a limit on the unabsorbed A1032A luminosity of 
-^1032 < 2 x lO^ergs" 1 if the O VI halo is co exten- 
sive with the X-ray halo. lEdgar fc Chevalier] ()1986D com- 
puted O VI luminosities for gas cooling from above 10 6 K 
(0.1 keV) as a function of M for isobaric and isochoric 
cases. The O VI limit then corresponds to a cooling rate 
of M < 2-3M yr- 1 . 

While this is substantially higher than the rate ex- 
pected in the 1-T model (0.4M© yr" 1 ), it is close to that 
expected in the 2-T model (~ 3M yr _1 ). If the true 
value is near the upper bound, then the fountain model 
is implicated, but a more sensitive observation could rule 
it out. It is also worth pointing out that the 2-T model 
leaves little room for a separate warm-hot halo compo- 
nent. 

4.1.2. Scale Height and Masses 

The coronal scale height may also distinguish between 
the 1-T and 2-T models. Specifically, we can measure 
the scale height in each component of the models assum- 
ing hydrostatic equilibrium (the "observed" scale height) 
and ask whether the value is consistent with our expec- 
tations based on the cooling and equilibrating times. If 
the observed scale height is inconsistent with these times, 
the model may be unphysical. 
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In hydrostatic equilibrium, n(z) — iiQe~ z / H , where H 
is the scale height. Since H = kT / ixm^g z is a function of 
temperature, we expect different scale heights for the 1-T 
and 2-T cases. If we assume the underlying gravitational 
potential is the same in both models, then we expect the 
0.2 keV 1-T model to have a scale height twice as large 
as the 0.1 keV component of the 2-T model. However, 
the 0.25 keV component of the 2-T model is responsible 
for most of the observed emission in the Chandra and 
XMM-Newton bandpass and would have a similar scale 
height to the 0.2 keV 1-T model. The surface brightness 
profile is not known to sufficient accuracy to distinguish 
between the models by itself. 

In the 1-T interpretation, the density scale height H 
is twice the measured X-r ay surface brightness s cale 
height zq (from this work or lBregman fc Pildislll994l ) . or 
H ~ 6 — 7kpc. For the same potential, and ignoring po- 
tential compli cating factors such as magnetic fields and 
cosmic rays (|Kalberla fc Kerdll998h . the 0.1 keV com- 
ponent of the 2-T model would be expected to have 
H ~ 3kpc while the 0.25 keV component would have 
H ~ 7kpc. Spectra extracted along and above the disk 
at various positions indicate that the 1-T model atmo- 
sphere is isothermal, so these values should be approxi- 
mately accurate. However, to verify these values we have 
measured the emission-weighted density as a function of 
z in a set of five strips (each 0.5 arcmin high) above and 
below the disk. We then obtain H from an exponential 
fit to n(z) = n e' z/H . 

To maximize the signal, we fit spectra extracted from 
both sides of the disk and both sides of the bulge simul- 
taneously. As noted above, the surface brightness profile 
(Figure 2]) indicates a different zq north and south of the 
bulge, but as the surface brightness is proportional to 
n 2 , small differences in H are exaggerated in zq. We fit 
the combined spectrum from each z with the 1-T and 
2-T models to obtain the emission measure in each case, 
freezing the metallicity in the 1-T model to its overall 
best-fit value and in the 2-T model to solar. We also 
freeze the absorbing column density at values estimated 
from the H I maps in the inner halo and at the Galactic 
value Njj — 6.5 x 10 20 cm~ 2 above lkpc. For the pro- 
jected depth of the extraction boxes we adopt 6kpc (the 
error incurred here goes as the square root of the depth, 
and our adopted value is based on the scale length of 
X-ray emission along f assuming azimuthal symmetry). 

For the 1-T model, we fi nd H = 5 ± 2 kpc , which is 
consistent with the 7 kpc of Brcgma n fc Pildisl (|1994l) as 
well as the 4kpc corona l scale height of the Milky Way 
(|Kalberla fc Kerd fl998l) . In the 2-T model we derive 
H = 6 ± 2 kpc in the 0.25 keV component and H = 3 ± 

1 kpc in the 0.1 keV component. The uncertainties are la 
errors based on fitting the emission- weighted densities, 
but do not include the systematic error from degeneracy 
in model parameters in the spectral fitting (which is more 
important for the 2-T model). For reference, we find 
n ~ 6x 10~ 3 cm~ 3 in the 1-T model, n ~ 4x 10~ 3 cm~ 3 
in the 0.1 keV component of the 2-T model, and no ~ 

2 x 10 -3 cm -3 in the 0.25 keV component. The scale 
heights correspond to X-ray masses (assuming the halo is 
cylindrical with a radius the size of the projected optical 
disk radius R ~ 6.5 k pc) of M ~ 3 x 10 8 M Q in the 1-T 
model (cf. 10 8 M Q in lBregman fe Pildisl[T99l and M ~ 



1 x 1O 8 M in the 0.1 keV component of the 2-T model. 
The 0.25 keV component also has M < 10 8 M Q . 

The agreement between these values and those derived 
from assuming an isothermal atmosphere with the same 
underlying g z above is encouraging, but we have ne- 
glected the cooling and equilibrating times. We now ask 
whether the "observed" scale heights in each model are 
consistent with these times. 

Hot plasma rising from the disk will reach equilibrium 
in about a dynamical time. The density and the temper- 
ature of the 0.1 keV component place it in t he subsonic 
or transonic regime (jHouck fe Br cgman 1990), so the the 
sound-crossing time is a reasonable proxy for the dynam- 
ical time (the free-fall time and the halo-filling time are 
similar). For kT = O.lkeV, c s ~ 0.08 kpc Myr -1 , and for 
kT = 0.25 keV, c s = O.lSkpcMyr" 1 . 

0.1 keV plasma injected at the sound speed takes about 
38Myr to rise 3 kpc, during which it will cool. The 
characteristic cooling time of a plasma with a density 
of the value near the disk is t coo i ~ 12Myr. We can 
also estimate i coo i from the mass of the O.lkeV com- 
ponent (10 8 M Q ) and a cooling rate of 3-MQyr" 1 from 
above to find ^33Myr. In either case, this is less than 
the crossing time. The clouds that form continue to 
rise to near the expected scale height before falling back. 
iHouck fc Bregmanl (|1990l ) show that under these circum- 
stances a dense layer builds up. It is possible that the 
surface brightness enhancement around the disk to the 
north (between 1 — 2 kpc off the disk) is this dense layer, 
but no corresponding layer is visible in the H I maps 
where the column density falls off monotonically. Also, 
the north side of the bulge only has about 30% more 
mass than the south in this region. This corresponds to 
less than 1% of the total coronal mass. 

The measured scale height in the O.lkeV component 
of the 2-T model appears to be inconsistent with what 
we would expect from gas actually rising at this tempera- 
ture. A substantial amount of 0.1 kev gas is also required 
above 3 kpc in the 2-T fit, but even if we put the fountain 
model aside it is unlikely that this gas could have cooled 
from the hotter medium. This is because the O.lkeV 
material will cool below X-ray visibility long before the 
0.25 keV gas cools to O.lkeV. We might also expect to 
see filamentar y structures in the X- ray maps if this were 
the case (e.g. iMcCourt et al.ll20lj ). but we do not see 
any (Section [2|). 

In contrast, the 0.2 keV plasma of the 1-T corona 
moves faster and cools more slowly and will equilibrate 
long before it cools. The 1-T model is also in bet- 
te r agreement with recent ga lactic f ountain mode l s such 
as lFraternali fc Binnevj (|2008ft and lMarasco et alj (|2012D 
that find an accreted corona is required to dissipate the 
angular momentum of material ejected from the disk. In 
these models, 10-20% of the total extraplanar gas mass 
must be accreted fr om the IGM. The ma ss of the H I 
halo is 1.2 x 1O 9 M (jOosterloo et al.ll2007f ). so the X-ray 
mass in the 1-T model is ^15-20% of the total (depend- 
ing on the scale height and assumed geometry), whereas 
the 0.25 keV component of the 2-T model comprises 4- 
8% of the total. The masses are unlikely to be much 
higher, so the 1-T model is apparently better able to 
provide an angular momentum sink. We note that the 
masses of other known components of the extraplanar 



18 



HODGES-KLUCK & BREGMAN 



TABLE 6 
Halo Component Masses 



Component 


Mass 


Depletion Rate 


Ref. 




(M Q ) 


(Moyr" 1 ) 




H 2 


fcwxlO 8 


~ 15 


1,2 


H I 


1.2 x 10 9 


20 - 40 


3 


HII 


fcwxlO 8 


- 15 


2 


Warm-hot 


1 - 2 x 10 s 


< 2 


4, this work 


Hot 


2 - 4 x 10 8 


< 2 


5, this work 



References. — (1) IGarcia-Burillo et all 11993) . (2) 
Howk & Savage! |200 C|K (3) lOosterloo et al.| 120071) , (4) 
Otte et all J2003I) . f 5") iBregman fc Pildisl J1994D 



Note. — Estimated masses in halo components discussed in 
Section 14.21 The depletion rate in the cooler components is as- 
sumed to be the total mass divided by the free-fall timcscalc 
(based on the scale height) and in the hotter components it is 
the cooling rate. There is no positive evidence for a separate 
warm-hot component, but we cannot rule it out. 

gas are smaller than the H I halo, with a combined to- 
tal in the cold molecular, wa rm neutral, and warm ion- 
ized phases of a fewx 10 8 M^ (| Garcia- Burillo et al]|l993 
IHowk fc Savag^[2000l) . 

4.1.3. Concluding Remarks 

If the hot halo can be described by a simple thermal 
model with a filling factor of unity, the cooling rates are 
on the cusp of ruling out the 2-T model, while the mea- 
sured scale height of each model is more consistent with 
the 1-T model. These arguments depend primarily on 
temperature and are less sensitive to the details of the 
physical model. For example, 0.2 keV plasma in a 1- 
T fountain model would still equilibrate before cooling. 
Thus, we conclude that the 1-T low metallicity model is 
better at explaining the existing data than the 2-T solar 
metallicity model, but these considerations do not prove 
an origin. The low metallicity in the 1-T model does not 
appear to be due to incorrect subtraction of the X-ray 
background or other systematic effects, but it is possible 
that the spectrum is undermodeled (cf. IStrickland et al.l 
120001 ). The 400 ks simulated XMM-Newton spectrum in 
FigurefTUl would rule out a 1-T model based on Fe-L lines. 
If the 2-T model were still disfavored, such results would 
point toward a more complex hot halo. It is also possible 
that the 2-T model represents a fountain+accretion sce- 
nario in which the hotter component (whose metallicity 
is poorly constrained, but somewhat above the O.IZq of 
the 1-T model) represents accreted gas that has equili- 
brated at a different scale height than the cooler fountain 
gas. The simulated XMM-Newton spectrum in Figure [TU] 
would strongly constrain the metallicity in the 2-T model 
if the 1-T model has already been ruled out. In this case, 
the accretion component would dominate farther from 
the galaxy. 

4.2. Origin of the Extraplanar Gas 

We now turn to the origin of the halo gas. In the clas- 
sic galactic fountain scenario, superbubbles break out of 
the disk, venting hot gas to the halo which cools and 
falls back as H I clouds. In this section, we will first 
demonstrate that this model cannot explain the extra- 
planar gas in NGC 891 if it is in the steady state. We 
will then consider alternatives, including cold accretion 



and prior outflows, before placing NGC 891 in the con- 
text of other massive late- type galaxies. We conclude by 
offering a speculative picture that NGC 891 and several 
other galaxies represent a period of enhanced star for- 
mation in the disk following a nuclear starburst that has 
since wound down. 

4.2.1. Galactic Fountains and the Steady State 

Considering the relevant timescales and disparate 
masses in the components of the extraplanar gas (Ta- 
ble [6]), it is evident that the atomic component could 
not have cooled from the hot corona in the steady state. 
This is because the cold (T < 10 4 K) gas free falls 
back onto the disk (in tg ~ 30 — 50Myr) at a rate of 
M ~ 30M Q yr _1 , whereas the total cooling rate from 
the X-ray and UV emitting gas is M < 3M Q yr" 1 . This 
rate could be as low as M < 0.5M Q yr _1 in the 1-T ac- 
cretion model if the only source of O VI luminosity is 
cooling coronal gas. In this case, t he instantaneous gas 
consumption rate by star formation (|Popescu et al.|[2004l ) 
is more than 4 times higher than the accretion rate. 

The halo components could exist in a steady state if 
the H I is not supplied by the hot gas, but rather rises 
into the halo via a (primar ily cold) galactic fountain (e.g. 
iFraternali fc Binnevl 1200(1 and variants). The viability 
of a cold fountain depends on its ability to expel enough 
gas and to keep it cold. In our view, if the fountain 
is powered by superbubble breakouts (the most natural 
way to launch material more than 1 kpc above the disk) , 
it fails to meet both criteria as we argue below. This 
rules out a halo formed by the steady-state action of a 
galactic fountain, hot or cold. 

First, it is worth laying out our superbubble breakout 
paradigm. Superbubbles form around OB associations 
as stellar winds and SNe evacuate a region in the sur- 
rounding ISM. Once a bubble has formed, i ts structure is 
thoug ht to be similar t o the p icture given inlWeaver et al. 
(1977), ISmith fc Coxl (|200l . and IMac Low fc McCrav 
(1988). Near the core, the hypersonic wind or super- 
nova blast wave freely expands. Farther out, the flow en- 
counters a stagnation shock and heats to X-ray emitting 
temperatures (10 6 < T < 10 s K). This hot medium fills 
most of the bubble. SNe that explode after the bubble 
has formed do not form remnants like isolated SNe, but 
instead convert their kinetic energy into thermal energy 
by heating the already hot interior via weak shocks. Ra- 
diative losses from the very hot material are small. The 
bubble is bounded by a shell of swept-up mass, which is 
thin early on but thickens in developed bubbles. The hot 
gas within evaporates mass from the inner layer of this 
shell via thermal conduction, and in a developed bubble 
the hot mass is dominated by evaporated gas. The shell 
itself is not significantly depleted by this process. 

Once a bubble's diamete r exceeds about thre e times 
the scale height of the disk (Mac Low et al.|[l989f ). it will 
break out into the halo. The breakout is characterized 
by accelerated upward shell expansion, which leads to 
Rayleigh- Taylor instabilities that disrupt the shell and 
allow the hot gas to vent. Howev er, less than 10% of th e 
shell mass erupts into the halo (|Mac Low et al.l fl989). 
and this rises into the halo as fragmented clouds. The hot 
gas evacuates as a pressure-driven outflow. This scenario 
finds some support in NGC 891: in a study of the emis- 
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sion fr om dust and diffuse Ha above the disk. lRossa et al.l 
(|2004f ) do not find large, organized "chimneys" . Rather, 
the evidence indicates that the bubbles vent through a 
disorganized medium. 

To determine whether the outflow is predominantly hot 
or cold requires estimating the masses and energies in- 
volved. We do not have direct observations in NGC 891 
(although dust filaments indicating bubbles and a some- 
what porous disk a re present), but we can estimate the 
cold mass from the iBoomsma et al.l (|2008j ) study of H I 
holes in NGC 6946, a nearby face-on spiral with a simi- 
lar mass and SFR as NGC 891. We also consult values 
derived from other face-on spirals that are not as appar- 
ently analogous to NGC 891. The hot mass is uncertain, 
but we can constrain it by considerations from nearby 
superbubbles that have been studied in the X-rays. The 
fiducial values we derive are given in Table [7j 

We assume that the cold mass ejected into the halo 
is 10% of th e mass of the swept - up shell at the time 
of breakout ()Mac Low et al.lll989i ). although it may be 
less. As will be clear later, it is difficult from an ener- 
getics perspective to expect much more. The breakout 
size of the hole depends on the scale heig ht of the disk, 
which in NGC 891 is about 200-400 pc (jXilouris et al.1 
Il998h . We adopt 300 pc as a fiducial value, meaning a 
bubble at breakout has a radius of about 450 pc, but we 
note that the sc ale height varies with d istance from the 
galactic center. IBoomsma et al.l (j2008| ) estimated shell 
masses in deep observations of H I holes in NGC 6946 
by measuring the ambient column density around each 
hole and the area of the hole. Because they measured 
Nh in regions including the swept-up shell itself, the re- 
sulting mass is an upper bound. They find an exponen- 
tial distribution of hole size with radius and an average 
missing mass per hole of M < 1O 7 M . This is consis- 
tent w ith the swept-up shell mass in the iMac Low et al.l 
( 1989) simulations of M ~ 3 x 10 6 M Q at breakout, and 
the lBoomsma et ail (|2008l ) value is biased towards higher 
mass by the regions used to measure Njj, the angular 
resolution (their synthesized beamsize is 390 pc, which is 
close to the breakout size), and hole identification algo- 
rithms (since M oc Apr 2 , larger hol e s hav e much more 



TABLE 7 
superbubble fiducial values at 
Breakout 



mass). For example. IBoomsma et al.1 (I200S ) iden tify 121 
holes in NGC 6946 while iBagetakos et al.1 poH find 56 
in the same galaxy. 

We can also compare these result s to the face-on 
galax ies M101 (lKamphuislll993l). M31 dBrinks fc Baiaial 
Il986jh and NGC 2403 dThilker et al.l 119981) where sim- 
ilar measurements have been made at high resolution 
and sensitivity. M101 has a similar mass and SFR 
to NGC 891 (~ 5Af Q yr _1 in IKuntz fc Snowded[20i0l ). 
whereas M31 has a similar mass but much lower SFR, 
and NGC 2403 is about ten times less massive but has 
a relatively high SFR for its mass (~ 0.6A/©yr _1 in 
iHeald et al.ll201lL see erratum) and a lagging H I halo, 
suggesting it is a lower mass analog to NGC 891. The 
results suggest that the amount of mass swept up for a 
shell of a gi ven size and age i s roughly the same across 
the galaxies. iKamphuisI (|1993|) find relatively large holes 
in M101 comparable to or larger than the supershells 
seen in NGC 6946, with similar masses and dyn amical 
ages. On the other hand, iBrinks fc Bajajal (|1986f ) gener- 
ally find smaller, younger holes in M31 than NGC 6946, 



Quantity 


Value 


SNe before breakout 


~100 


Energy budget 


~ 10 53 erg 


Shell mass 


3 x 1O 6 M 


Disk scale height 


300 pc 


Radius at breakout 


450 pc 


Breakout cold mass 


3 x 1O 5 M 


Breakout hot mass 


5 x 1O 4 M 


n e 


0.01 cm" 3 


T e 


4 x 10 6 K 



Note. — Fiducial values used in 
Section 14.21 We consider these rea- 
sonable estimates of the true values, 
which arc hard to determine. See text 
for details. 



with a typical missing mass of only 1O 5 M per hole for 
holes with radii of abo ut 100 pc and kinem atic ages of 2- 
30Myr. In NGC 2403. [Thilker et al.l ([T998T ) identify both 
kinematically young (expanding) shells and older stalled 
shells with a missing mass range of 10 6 — 10 7 M Q per hole 
and ages between 1-200 Myr. The range of bubble sizes, 
ages, and masses generally supports the IMac Low et al.l 
(1989) picture and indicates that bubbles large enough to 
break out of the disk in NGC 891 would have masses be- 
tween 10 6 — 10 7 M Q . Based on the distribution of the hole 
sizes and masses in the four face-on galaxies we consider, 
we adopt a fiducial missing mass of a r = 450 pc hole in 
NGC 891 to be M = 3 X 10 6 M o . The cold mass expelled 
upwards at breakout would the n be M < 3 x 1 5 M ( ti 
per hole (|Mac Low et al.l 11989ft . IMac Low etaLl (|1989f i 
find that the ejected shell material has upward velocities 
between 50-100 km s" 1 , so the kinetic energy required is 
7 x 10 51 erg< E < 3 x 10 52 erg. This is comparable to 
the total kinetic energy in the rest of the shell, but small 
compared to the thermal energy stored in the hot gas 
that drives the outflow (see below). 

The hot mass in bubbles is quite uncertain, but we 
estimate the mass by assuming a spherical bubble just 
prior to breakout whose size is 3 times the scale height 
of the disk in NGC 891. A bubble filled with hot gas 
has a mass M/f ~ ^irr 3 /j,mHn e , where / is the fill- 
ing factor. Although we expect the bubble to be filled 
with hot material (the sound speed of the hot gas is 
much faster than the bubble expansion rate), much of the 
(soft) X-ray emission comes from the inner e dge of the 
shell where supernova bl ast waves impact it (Ch u et aLl 
119951 : Uaskot et all 1201 lh . The rest of the bubble is 
presumably filled with even hotter, more tenuous gas. 
Thus, the emission measure of the hot component that 
dominates the mass is n%fV. Electron densities mea- 
sured in X-ray observations of nearby superbubbles range 
from n e \/J ~ 1 — 100 x 1 0~ 3 cm" 3 for bubbles between 
100 - 600 pc across (e.g ICash et all 119801: IPoints et al.l 
120001: iDumie et al.ll200lt ISmith fc Wand 120041) . For the 
breakout bubble radius of 450 pc, the mass is M/f ~ 
5 x 10 5 MQ/ij^ pc n_3, where ft-kpc is the disk scale height 
in kpc and n_3 is in units of 10~ 3 cm~ 3 . For fidu- 
cial values of h = 0.3 kpc and n — 0.01cm -3 , we find 
M//-1O 5 M . 
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There is considerable uncertainty regarding the filling 
factor of the X-ray emitting material in nearby super- 
bubbles, and we do not have direct information from 
NGC891. However, we suggest that a reasonable lower 
bound for the mass of hot material in a bubble at break- 
out is M ~ 1O 4 M , which we justify presently. Over 100 
SNe are thought to be required for a bubble to break 
out of the disk of a galaxy like the M ilky Way (e.g. 
iMac Low et~aT1ll989t iMelioli et aTll2008D . and by exten- 
sion NGC891. This will release about 10 53 erg (stellar 
winds will also contribute to the kinetic luminosity). Of 
this, a few x 10 52 erg is contained in the average shell de- 
scribed in Booms ma et al.l (|2008f) . The majority of the 
energy is either radiated or goes into heating gas evapo- 
rated from th e inner wall of the shell to X-ray emitting 
temperatures (jWeaver et al.lll977l ). Nearby superbubbles 
have interior temperatures ranging from 1 — 8 x 10 6 K 
(e.g. [Cash et alJll980HDunne et al.ll200H ISmith fc Wanal 
l200l . so as a fiducial temperature we take T = 4 x 10 6 K. 
The specific energy of material at this temperature is 
e = 8 x 10 14 ergg-\ or e = 1.6 x lO^ergM" 1 . If about 
30% of the superbubble energy budget goes into heating 
evaporated gas, about 2 x 1O 4 M of gas could be heated 
to the fiducial temperature. The actual value may be 
higher depending on radiative losses. 

Extrapolations from observed superbubbles also sup- 
port M ~ 10 4 M Q of hot gas as a lower limit in large bub- 
bles. If the hot m asses inferred in smaller bubbles (e.g. 
Dun ne et al.ll200l are scaled up to bubbles of r = 450 pc, 
the masses are M ~ 1 - 10 x 1O 4 M , and lYukita et al.l 
(2010) find in spectral fits to H II regions in NGC 2403 
that the ratio of hot gas to swept-up mass is 6-10% (for 
our fiducial shell mass of M = 3 x 10 6 M oc ; ot , this is 
M - 18 - 30 x 1O 4 M of hot material, assuming this 
ratio holds throughout the lifetime of a massive bubble) . 

These considerations suggest a hot mass somewhere 
between 1 — 10 x 10 4 M Q , so for the sake of argument 
we take as a fiducial value the hot mass in a r = 450 pc 
bubble to be M = 5 x 1O 4 M . This would give a filling 
factor / = 0.5 in the bubble described above, but we 
emphasize that the hot mass and filling factor of the soft 
X-ray emission is poorly known in many bubbles. 

Thus, at breakout our fiducial bubble (Table [7]) 
launches 3 x 1O 5 M of cold shell material towards the 
halo and contains 5 x 1O 4 M of hot (T = 4 x 10 6 K) ma- 
terial whose thermal energy (E = 8 x 10 52 erg) is similar 
to or greater than the kinetic energy in the entire shell. 
The erupting shell material is therefore driven rapidly 
upwards even though it dominates the mass. Since the 
hot gas flows faster than the ejected shell fragments, it 
can evaporate the neutral material just as in the con- 
fined bubble. At our fiducial temperature of 4 x 10 6 K, 
the sound speed is c s ps 150 km s -1 , whereas the shell 
ejecta travel at 50 — 100 km s -1 in the simulations of 
IMac Low et al.l (|1989l ). Thus, the cool parcels can be 
considered as islands in a stream. The thermal energy 
of the hot gas is sufficient to raise about 2 x 1O 5 M of 
cold material to 10 6 K, so it is energetically feasible for 
the hot gas to become the dominant component (there is 
more than enough energy to merely evaporate the cold 
material by heating it to, say, 10 5 K). 

Whether the hot component becomes dominant is un- 



clear, since evaporation occurs via thermal conduction in 
a boundary layer around the fragments and most of the 
hot material would not participate. On the other hand, 
the fragments are rising into the hot bath of the extant 
halo. A significant unknown is the effect of magnetic 
fields, which could suppress conduction relative to the 
iSpitzerl (|1962f) value if the fragments are wrapped in field 
lines, which are seen to be expelled along with the gas 
(|Healdll2012l) in at least one superbubble in NGC 6946. 
However, t urbulence could mitigate this effect (cf. galaxy 
clusters in iNaravan & Mcd vedevl 120011 ) . In the bubble 
described by iHealdl (|2012h . both the large-scale fields 
around the bubble and turbulent fields entrained in the 
hot plasma are uplifted. The speed of the hot gas flow is 
also important, since gas that is too fast shocks the cooler 
entra ined gas instead of evaporating it ({Strickland et al.l 
12000( 1. Pure hydrodyna mic models suggest that virtually 
all of the gas rises hot (jMelioli et al.1 120081 ) . but the de- 
tection of polycyclic aromatic hydroc arbons (PAHs) to 
a scale height of -500 pc in NGC 891 (jRand et al.ll201lD 
may indicate that a substantial amount of material rises 
cold. The difference between PAH features in the disk 
and halo suggests that the cool fountain gas does expe- 
rience shocks as it rises, but that these do not heat the 
gas to X-ray emitting temperatures (it is also possible 
the grain destruction time exceeds the flow time). 

The proportion of cold fragments to hot gas will de- 
cline with distance above the disk since some of the cold 
material that is not evaporated will fall back below 1 kpc 
(above wh ich it would unambigu ously be considered halo 
material). IMac Low et al.l ([1989) find that the fragments 
have upward velocities between 50 — lOOkrns -1 , and on 
the lower end these fragments will reach an apex below 
1 kpc. 

Considering these factors, we suppose that less than 
half of the fountain material that makes it into the halo 
is cold. This amounts to a fewxlO 5 M of cold material 
injected into the halo per hole and a comparable mass of 
hot gas. Thus, a galactic fountain powered by supernova 
feedback seems unlikely to be truly "cold" . 

More importantly, it is evident from these values that 
the proportion of cold and hot material in the fountain 
is a moot point in NGC 891 because the fountain envi- 
sioned above does not expel enough mass to sustain the 
halo regardless of the temperature. This can also be seen 
in NGC 6946, whose H I holes we used to derive our val- 
ues. For our fiducial values, the total amount of gas that 
rises into the halo in NGC 6 946 is no more than — 15% 
of the total extraplanar H I (|Boomsma e t aL 2008]) and 
coronal gas (jSchlegel et al.l 120031 ) . Even with optimistic 
assumptions about the true hole size and amount of ki- 
netic energy involved, it is difficult to envision increasing 
the outflow mass by an order of magnitude without eject- 
ing the entire swept-up shell mass of each contributing 
hole, which is impossible for isolated holes. The situa- 
tion is even worse in NGC 891, which has a similar SFR 
to NGC 6946 (and therefore presumably about the same 
number of holes) but four times more extraplanar H I. 

This problem is circumvented if the extraplanar gas 
spends more time in the halo than the lifetime of the 
progenitor bubble. The lifetime of the H I holes is lim- 
ited by the slowing expansion velocity (once it drops be- 
low the characteristic velocity in the surrounding gas, 
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the bubble breaks up) and interactions with neighboring 
bubbles (converging flows from two expanding shells can 
promote star formation which w ill disrupt that part o f 
the shell), and is perhaps 80Myr (jBoomsma et al.l f2008). 
Hot gas may reside in the halo longer than this, but we 
have already established that the H I halo cannot have 
cooled from the X-ray emitting gas if it is in the steady 
state. The time spent in the halo by the cold gas is about 
twice the free-fall time, or 30 — 80Myr for the ejection 
speeds in IMac Low et all (| 19891) . Since the expulsion of 
cold mass from the bubble is a one-time event occurring 
at breakout, we expect the cold material launched into 
the halo to return to the disk within a bubble lifetime 
(although it will not refill the hole since the galaxy ro- 
tates and the fount ain material has angular momentum; 
IMelioli et al. I [2008V Additional hot material may leak 
out of the hole following the breakout, but for an isolated 
hole this will be subsequent stellar wind and SNe mate- 
rial. The masses of these components is quite small; even 
if we assume that the hot gas does not escape the blown- 
out bubble before evaporating additional material from 
the shell, the evaporated component described above was 
a small fraction of the total swept-up shell mass. 

Thus, neither a classic nor a cold fountain seem ca- 
pable of producing the extraplanar H I in NGC 891 in 
the steady state at the current SFR. Although our ar- 
gument for NGC 891 is largely circumstantial, we note 
that NGC 6946 has a similar mass, SFR, and is disk- 
dominated like NGC 891 with a comparable disk scale 
height. Since stellar feedback must be tied closely to the 
SFR, there is no reason to expect a much higher hole 
rate in NGC 891. 

It is not clear if this imbalance occurs more generally. 
There are not many non-starburst galaxies with the X- 
ray, UV, Ha, H I, and radio continuum data required 
to establish an H I mass and cooling rate (starbursts 
are not expected to be in equilibrium). This is largely 
due to the paucity of X-ray and O VI m easurements 
aroun d non-starburst massive spirals (see I Tyler et al.l 
I2004h . However, among galaxies with relatively com- 
plete coverage, there are others that appear to fall 
into this category. As we have established, the H I 
halo of NGC 6946 cannot have been ejected by the ob- 
served holes unless the expelled mass rate is an order 
of magnitude larger, but it also cannot have cooled 
from the hot corona (which has a bolometric Lx < 
lO^ergs" 1 : iSchlegel et all 12003ft . M101 likewise has a 
similar mass and SFR with a bolometric halo X-ray lu- 
minosity of L x < lO^ergs" 1 (|Kuntz & Snowdenll2010D 
with H I holes detected throughout the disk (Kamphuis 
1993). The halo X -ray cooling rates implied are M < 
O.O5M yr _1 . NGC 2403 has the same issue: the X-ray 
emitting mass c onceivably in the ha lo is only 10% of the 
extraplanar H I (|Yukita et al.H2010H . 

Unfortunately, the picture is not clear in normal 
edge-on spirals. NGC 891 has the highest X-ray 
surface brightness of a non-starburst edge-on galaxy 
(jBregman fc Houckl 119971) . so much of the X-ray work 

on edge-on spirals has focused o n starbursts . A 

few exceptions i nclud e NGC 4631 (IWang et al.l 119951: 
Strickland et~aH [200l and NGC 4565 dVogler et al.1 
1995ft. but both are in teracting with companions, and 
Strickland et al.l ()2004D argue that NGC 4631 is experi- 



encing a weak starburst. NGC 4565 has a similar mass 
to NGC 891 but a much smaller SFR (~ OJAfayr" 1 in 
iHeald et all 1201 it see erratum) and evi dence only for a 
weak H I halo dZschaechner et all 1201 lh . NGC 4631 has 
a large halo, but IRandl ()1994f ) argue based on kinemat- 
ics that most of the H I far from the disk is due in part 
to the not perfectly edge-on inclination (85°), tidal in- 
teraction, a nd a large scale height in the outer part of 
the galaxy. IStrickland et al.l (|2004D included two non- 
starburst edge-on galaxies in their Chandra sample aside 
from NGC 891 (NGC 6503 and NGC 4244), but both are 
substantially less massive than NGC 891 and si gnificant 
X-ray emission was detected in neither. Finally. iLi et al.l 
(12007ft detected diffuse X -ray emission in M104 with a 
bolometric X-ray luminosity of a fewxl0 39 ergs -1 , but 
M104 is also at least s evera l times more massive than 
NGC 891 and ILi et all (12007ft argue that the halo is X- 
ray under-luminous while iBaiaia et all (|1984l ) find a rel- 
atively small amount of H I in the galaxy. 

If the halo is not in the steady state, there are two 
basic possibilities for the origin of the extraplanar gas: 
cold accretion and prior outflows. 

4.2.2. Cold Accretion 

Cold accretion (including the recycling of gas ripped 
out of a disk by galactic interaction) can occur via 
the infall of intergalactic H I clouds or through tidal 
stripping and cannibalization of small satellite galaxies. 
The former is disfavo red by the paucity of intergalactic 
H I (see discussion in lSancisi et all 12008') and the small 
amount of counter-rot ating extraplanar H I in NGC 891 
(jOosterloo et al.ll2007f ). Recent simulations also indicate 
that accret ion of fresh material primarily occurs in th e 
hot phase (jQppenheimer et al.ll2010l : Uoung et al.| [2012) . 
While it is possible that a large reservoir of photoionized 
hydrogen resides outside the visible H I halo, this would 
requir e an extremely extended halo. The Oost erloo et all 
(2007) contour map extends down to 10 19 cm -2 and 
monotonically declines from the center. The scale height 
suggests that the column drops below 10 18 cm~ 2 within 
several kpc of the last contour. At this point, we would 
expect the H I to be photoionized. Although the column 
density likely continues to drop beyond this, if we assume 
that the space around NGC 891 is filled with a uniform 
column of 10 18 cm -2 (10 4 M© kpc -2 ), it would take about 
10 5 kpc 2 to make up 10 9 M Q of ionized hydrogen. If we 
further assume that the photoionized halo is circular in 
projection and centered on NGC 891, the radius required 
is about 200 kpc. Replenishing the visible H I halo from 
such a reservoir cannot explain the observed mass, and it 
is not clear that such a halo would rotate with the galaxy 
to such distances. 

A more plausible scenario is that the halo is made up of 
tidally stripped gas from a satellite galaxy and gas ripped 
out of the disk of NGC 891 by the same interaction. How- 
ever, there is no direct evidence of an interaction such as 
tidal streams, stellar shells, or a warped disk that would 
be particularly visible in an edge-on galaxy. The absence 
of these features and the kinematics of the extraplanar 
H I indicate that the time elapsed since the posited in- 
teraction is a few times the galactic rotation timescale of 
a fewx 10 8 yr. In this case, however, we would expect gas 
within a few tens of kpc to have fallen back onto the disk. 
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There are a few candidates for recent interactions. First, 
the companion galaxy UGC 1807 is a small but gas-rich 
galaxy with about 10% of the mass of NGC 891 and is 
suffici ently far away to acco unt for the l ack of tidal fea- 
tures (jOosterloo et al.l 12007). However. lOosterloo et al.1 
(2007) find that the amount of accretion that can be ex- 
plained by this interaction is only M > 0.1M Q yr -1 . 

Two other sig natures of recen t or o ngoing accretion 
are described in iMouhcine et al.1 (|2010T l , namely an ex- 
tended (lOkpc above the disk), flat (disk-like) enve- 
lope of halo stars and a giant rosette pattern of tidal 
streams extending out to 30kpc from NGC 891. We 
do not believe the former can account for much recent 
cold accretion because such extended "disks" form on 
long timescales of a few Gyr or more when dwarf galax- 
ies are cannibal i zed (|Penarrubia et al.l[2006h . Moreover, 
IMouhcine etHI pnof T also find a similar extended halo 
around NGC 2683, a highly inclined (80°) galaxy with a 
similar total ma ss to NGC 891 but no substantial H I halo 
(Kennedy 2009). Such envelopes may be common, but 
evidently need not correla te with recent cold acc retion. 
Regarding the rosette tails. IMouhcin e et al.1 (|2010t ) argue 
that they likely originate from a single accretion event 
based on the similar morphology to disrupted dwarf 
galaxies in A-body simulations. If the galaxy were gas- 
rich, it could provide a substantial col d accretion rate 
over a relatively long period of time, but IMouhcine et"afl 
(2010) argue based on stellar ages (no stars younger than 
a few Gyr were seen in the streams) that the compan- 
ion had little H I compared to UGC 1807. Although 
the instantaneous cold accretion rate could be quite high 
upon absorption of a particular stream, we would expect 
an average accretion rate from a gas-poor dwarf with 
a pre-merger gas mass Mm ~ 10 8 M Q to be no more 
than M ~ 0.1 — 0.2AfQyr _1 , considering that the life- 
time of the stellar stream s is between 0.5 — 1 Gyr (e.g. 
iBinnev fc Tremaind l2008f ) . This is similar to the im- 
puted contribution fro m UGC 1807 of M > O.lAf© yr" 1 
(jOosterloo et al.ll2007l ). 

4.2.3. Prior Outflow 

The remaining possibility we consider is a prior hot 
outflow that cooled to form the H I halo. There are two 
basic models, including a vigorous fountain (powered by 
higher star formation distributed throughout the disk) 
or a nuclear wind. Plausible cooling mechanisms that 
would produce the observed halo can be found for either 
case, but as winds a re ubiquitous in starburst galaxies 
(jVeilleux et al.ll2005h , it may be that a powerful fountain 
is always accompanied by a nuclear wind. 

Typical outflow rates in winds a ppear to be ~ 2 — 3 
times the instantaneous SFR (e.g. iBouche et all [2012). 
although it is unclear how much mass they carry because 
the mass carried in each compon ent is hard to measure 
(see lStrickland fc Heckm an 2009). If material flowed out 
of the disk at 3OM yr" 1 , this would imply a SFR of 10- 
15Mq yr -1 during the active phase, which is plausible for 
massive starbursts. 

As the putative outflow must have occurred long 
enough ago for the H I to have cooled from the hot phase, 
a variety of parameters are possible. The minimum age 
of the outflow can be constrained by the cooling time of 
the gas and (in the case of the wind) the fallback time 



from large distances. In either case, ages of 50Myr or 
more are expected. The abs ence of wind signature s (cf. 
edge-on starburst galaxies in iStrickland e t al. 2004]) and 
the duration of starburst episodes place lower limits con- 
sistent with this value. Thus, it is neither possible to 
rule this out nor to make verifiable predictions, so what 
follows is speculation. 

NGC 891 belongs to a class of massive, disk-dominated, 
late-type galaxies with relatively high disk-distributed 
SFR that are not classified as starbursts. Nearby exam- 
ples include NGC 6946, NGC 2403, M10 1, NGC 4631 (al- 
though it is classified as a starburst bv IStrickland et al 



2004), and NGC 5775 (classified as a starburst in 



Li et al 



20081) . Of these, NGC 4631 and NGC 5775 show some ev- 



idence for a wind-like nuclear outflow in the X-rays, with 
the case being stronger in NGC 4631. As one might ex- 
pect, such galaxies tend to have more prominent halos 
than other massive spirals with lower SFR (e.g. iRandl 
Il996t lDahlemlfl9"97j) . with an active fountain implicated 
as a result. The existence of H I halos that appear too 
massive to have cooled in the steady state in NGC 891, 
NGC 6946, M101, and NGC 2403 (more precisely, the 
X-ray cooling rate and expected cold outflow appears 
unable to replace gas lost via free-fall) suggests that 
NGC 891 may simply be an extreme case of a more fre- 
quent phenomenon, and that prior outflows may also be 
necessary in these galaxies. The case of NGC 2403 is 
especially interesting in this vein because a Chandra ob- 
servation shows X-ray emission that is m ore centrally 
conce ntrated than the star-forming regions (lYukita et al.1 
120101) . Since NGC 2403 has a relatively high SFR (~ 
lMoyr^ 1 ) for its mass (M T ~ 10 10 M Q ) and the X- 
ray emission (if produced by star formation) represents 
a significant amount of mass, NGC 2403 may have had a 
nuclear starburst that has propagated out wards or shut 
off w hile the disk star formation continued (jYukita et al.l 
120101 ). Another hint of a prior outflow is found in 
M104, which h as a curiously low SFR for its mass, where 
ILi et all (pOOl argue that a nuclear wind (in their model 
powered by Type la SNe) has removed much of the ma- 
terial in the inner region of the galaxy, an argument 
supported by the absence of H I in the same region 
(Bai aia et alJll984h . 

We therefore speculate that the high SFR is the re- 
mainder of a starburst phase that ejected a substantial 
amount of material or else fresh star formation fueled by 
the return of this material to the disk. This scenario is 
attractive because it explains the evidence that points 
towards a vigorous ongoing disk-halo interaction (sum- 
marized in Section [T] for NGC 891), the apparent discon- 
nect between the H I and X-ray halos (and perhaps the 
low metallicity of the latter, if it is accreted), and why 
starbursts generally do not have extended hydrostatic 
X-ray halos or prominent extraplanar H I (although the 
latter may be due to photoionization of extant hydro- 
gen). This appears to be true in NGC 253, NGC 3628, 
NGC 3079, and the less massive NGC 4945, NGC 1482, 
and M82, where outflow m orphology is seen in the X- 
rays ((Strickland et al.ll2004D an d other bands , but only 
modest signs of extraplanar H I (Dahlcm 1997). The out- 
flow rates required to explain the H I halos are easily met 
in starbursts. On the other hand, if the H I halos develop 
as hot gas cools, we would expect to see X-ray luminous 
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post-starburst galaxies. This seems likely even if much 
of the H I halo is formed by recombination of extant 
photoionized gas in the wake of the starburst (cf. the 
starburst galaxy NGC 253 where the H I dis k is smaller 
than the optical disk; Boomsm a et al.ll20Q5D . Although 
the sample of nearby massive spirals with X-ray obser- 
vations is small, NGC 891 is the most X-ray luminous of 
the nearby non-starburst spiral galaxies. 

In principle, this picture is consistent with either the 
1-T or 2-T coronal models depending on the type of out- 
flow. It is possible that the observed corona is dominated 
by accreted material that has fallen in during and since 
the formation of the H I halo. Such a situation could arise 
in an outflow scenario if a distributed starburst (i.e., non- 
nuclear) "fountain" expelled a large amount of material 
(~ 1O 9 M ) in a relatively short period of time (a few 
tens of Myr) . If the material was expelled via superbub- 
ble breakouts at a constant rate and entered the halo at 
several million degrees, the initial material would equili- 
brate before the end of the starburst and begin to cool. 
Although its cooling time would be quite high at first, 
the scale height of the corona would be determined by 
the temperature, so additional fountain action would in- 
crease the coronal density (and luminosity), thereby de- 
creasing the cooling time. As more material accumulated 
in the halo, the cooling time would eventually shrink un- 
til the cooling rate was comparable to the fallback time of 
cold material, and the coronal gas would rapidly cool to 
form a large H I halo. If there were steady accretion from 
the IGM throughout this period, the hotter material fur- 
ther out would survive the rapid cooling phase, leaving 
the halo possibly dominated by accreted gas. Whether 
this would actually happen depends on the fountain ejec- 
tion rate, the temperature, and the details of the cool- 
ing. For example, H I clouds may primarily form via 
the thermal instability, and the putative "cooling flow" 
would form the most cold gas nearest the disk, so the 
H I halo might not extend more than a few kpc from the 
disk. More importantly, it is not clear that such powerful 
distributed "starbursts" occur. On the other hand, it is 
easier to explain such an outflow in the 2-T model of the 
halo (if the 0.25 keV component has near solar metallic- 
ity and the 0.1 keV component represents the present-day 
fountain). 

Overall, the hypothesis of a previous outflow is severely 
hampered by the dearth of meaningful constraints, but 
the present-day fountain and cold accretion do not seem 
sufficient to explain the massive H I halo in a steady 
state. 

4.2.4. Concluding Remarks 

We have argued that the massive H I halo in NGC 891 
(and, we speculate, in a handful of other nearby galax- 
ies with sufficient data) cannot be supported by the in- 
stantaneous SFR or cooling from the X-ray halo. How- 
ever, a rigorous examination of the fountain models and 
halo masses must be made in order to confirm this pic- 
ture. For example, it is possible that the interaction be- 
tween multiple brea kouts within a small patch of the disk 
(jMelioli et al.ll2009l ) might conspire to eject much more 
mass, or that new bubbles forming below the falling back 
fragments of old ones can lift much more cold material 
than the breakout of a bubble from a pristine disk. Our 
arguments are based on the latter scenario as it appears 



in maps of H I holes that the overlap fraction is small, 
but this may be a bias of the hole detection algorithm. 
An additional and important unknown is the amount of 
warm-hot material in the halos of various galaxies, al- 
though we would expect a separate warm-hot component 
to equilibrate with a small scale height. In NGC 891, 
at least, such a component could not explain the sub- 
stantial cxtraplanar H I a few kpc away from the disk. 
Finally, we emphasize that any galaxy with a low lumi- 
nosity (Lx < 10 39 crgs _1 ) X-ray halo (including many 
nearby massive spirals) must have a low cooling rate of 
M < lM©yr -1 , with M < O.lM yr _1 for a tempera- 
ture like that of the Galactic halo. 

Even if t he fountain model ca n be reconciled with 
NGC 6946 (jBoomsma et al.lf2008h . NGC 891 has a H I 
halo four times more massive with a similar instanta- 
neous SFR. The mass loading would then need to be a 
factor of 4 — 8 in order to balance the fallback rate of the 
extraplanar H I. If the H I within 1 kpc of the disk (but 
outside a few scale heights) is also participating in the 
fountain, the mass loading factor must exceed about 10 
(this may be material that was not ejected with sufficient 
velocity to exceed 1 kpc) . It is difficult to envision any 
fountain scenario in which this could be the case. 

Since cold accretion does not appear to be a viable 
source of H I, we speculate that a prior outflow is the best 
explanation. Unfortunately, there are few constraints on 
the characteristics of such an outflow, so it is an unsatis- 
factory, if plausible, explanation at present. However, if 
the present activity in the disk is insufficient to produce 
the neutral halo, it is the best candidate. Finally, if it 
is conclusively demonstrated that the X-ray emitting gas 
is indeed low metallicity, a prior outflow scheme would 
decouple it from the observed H I halo (that presumably 
has the same metallicity of the disk). 

5. SUMMARY & CONCLUSIONS 

We have used deep X-ray exposures of the nearby, 
edge-on spiral galaxy NGC 891 to determine the prop- 
erties of its hot halo. Our main findings are summarized 
here: 

• The hot halo is fit well by a thermal model with 
subsolar metallicity (Z ~ 0.1Z@) to high confi- 
dence. Fits to the halo emission near the disk and 
farther out find nearly identical parameters, indi- 
cating a homogenous halo. 

• The simplest explanation for a low- metallicity X- 
ray halo is hot accretion from the IGM. 

• The best-fit values agree between the MOS, pn, and 
Chandra ACIS-S detectors. An examination of the 
X-ray background in the region indicates that it is 
properly subtracted from the spectra, and known 
systematic effects from model choices and fitting 
method do not substantially change our results. 

• A 2-T model with solar metallicity is also a good 
fit to the halo spectrum when the emission measure 
is dominated by an absorbed, cool, kT ~ 0.1 keV 
component. Most of the emission, however, comes 
from a weaker kT ~ 0.25 keV component. The 
metallicity in this model supports a disk origin for 
the coronal material. 
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• When considering the X-ray cooling rate and O VI 
limits on the cooling mass as well as the density 
scale heights in the 1-T and 2-T models, the 1-T 
model is favored. 

• The aforementioned cooling rates are not high 
enough to offset the losses of the massive H I halo 
due to fallback of cold material onto the disk. This 
suggests that the extraplanar gas cannot be in equi- 
librium. 

• An examination of single superbubble breakouts in- 
dicates that a galactic fountain cannot lift enough 
material into the halo to account for the massive 
H I halo. This corroborates the cooling rate prob- 
lem. 

• There is no plausible cold accretion mechanism 
to supply so much gas so rapidly to NGC891, 
so we speculate that a prior outflow has taken 
place. There is circumstantial evidence for some 
kind of past nuclear starburst in a few other mas- 
sive, late-type galaxies with high star formation 
spread throughout the disk. However, the parame- 
ters of such an outflow cannot be constrained from 
the present data. 

1-T and 2-T models of hot halos have been criticized as 
undermodeling the tr ue spectrum, which may be consid- 
erably more complex ((Strickland fc Steven s 2000) . How- 
ever, there is reason to believe (at least in non-starburst 
galaxies) that the halo really is simple. Although we 
have presented arguments in favor of the reality of the 
low Z in the 1-T model and that the 1-T model better 
explains the coronal morphology, a definitive test could 
be made with a deep (400 ks total) exposure that would 
allow line ratios to be measured even at the low spec- 
tral resolution of the EPIC CCDs. We believe that we 
have demonstrated that the background is well enough 
understood and the temperature well enough established 
to make this test feasible. 



In addition, such an exposure would allow a definitive 
study of the smoothness of the halo. Coronal material 
is frequently described as "filamentary" in the literature, 
but in many cases the appearance of filaments may be 
biased by adaptive smoothing or related to PSF wings. 
In our XMM-Newton images of NGC 891, some apparent 
filamentary structure is revealed to be slightly brighter 
regions of the halo connected by dim point sources (seen 
in the Chandra data) or arcs that clearly trace the edge 
of a subtracted point source mask. Since the presence of 
this filamentary structure is an important test of state- 
of-the-a rt cooling models for accreted or hydros tatic ha- 
los fe.g. lMcCourt et al.|[20TTI: iHobbs et alj|2012j) . a deep 
search for filaments that are not obviously related to out- 
bursts (as in starburst galaxies) would be a welcome con- 
tribution. 

Progress on the steady-state front requires a more so- 
phisticated treatment of the galactic fountain mecha- 
nism, allowing for the action of interacting superbub- 
ble breakouts and superbubbles that form beneath de- 
bris falling back from prior breakouts. If it still ap- 
pears that the instantaneous SFR in NGC 891 and its 
high SFR counterparts cannot explain the amount of ex- 
traplanar gas, a serious search for a sequence of post- 
starburst galaxies with luminous X-ray halos is war- 
ranted. Improved H I im aging from the HALOGAS 
project ijHeald et al.l 120111) will refine the estimates of 
H I residing above 1 kpc in nearby galaxies, but unfor- 
tu nately X-ray imag ing r emains sparse. Sti l l, the work 
bv lTvler et all (pOOl and lKuntz fc Snowdenl (pOlOt) indi- 
cates that most galaxies will not have X-ray halos nearly 
as luminous as NGC 891. This in turn places cooling rate 
limits of about lM Q yr _1 on the halos of these galax- 
ies, assuming a single-temperature plasma in hydrostatic 
equilibrium. 

The authors thank the referee for his or her careful at- 
tention to this work and for helpful suggestions that im- 
proved the discussion of the H I halo. EHK gratefully ac- 
knowledges support from NASA ADAP grant #061951. 
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